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cn ' Abstract 

> 

il^ ' We derive super symmetric quantum mechanics of n BPS objects with 3n po- 

[-^ . sition degrees of freedom and An fermionic partners with 50(4) R-symmetry. 

^^ \ The potential terms, essential and sufficient for the index problem for non- 

C"~~- ' threshold BPS states, are universal, and 2(n — 1) dimensional classical moduli 

spaces Mn emerge from zero locus of the potential energy. We emphasize 
that there is no natural reduction of the quantum mechanics to A^^, contrary 
to the conventional wisdom. Nevertheless, via an index-preserving deforma- 
tion that breaks supersymmetry partially, we derive a Dirac index on A4„ as 
the fundamental state counting quantity. This rigorously fills a missing link 
C^ \ in the "Coulomb phase" wall-crossing formula in literature. We then impose 

Bose/Fermi statistics of identical centers, and derive the general wall-crossing 
formula, applicable to both BPS black holes and BPS dyons. Also explained 
dynamically is how the rational invariant ~ Q{l3)/p'^, appearing repeatedly in 
wall-crossing formulae, can be understood as the universal multiplicative factor 
due to p identical, coincident, yet unbound, BPS particles of charge (3. Along 
the way, we also clarify relationships between field theory state countings and 
quantum mechanical indices. 
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1 Introduction 

Wall-crossing refers to a phenomenon where BPS spectrum [1, 2] changes abruptly as 
either a parameter or a vacuum of a supersymmetric field theory is smoothly changed. 
In the four dimensional setting, it was first discovered in the context oi SU{2) Seiberg- 
Witten theory [3, 4], where the asymptotic spectrum consists of massive vector mesons 
and an infinite tower of dyons with unit magnetic charge while the spectrum at the 
center of the vacuum moduli space is composed of a monopole and a dyon [5] . This 
apparent "decay" or "disappearance" of BPS states occurs across a wall of marginal 
stability in the moduli space, and this phenomenon turned out to be generic in all of 
N >2 D = 4 theories. 

The mechanism of the disappearance was clarified a few years later. It was found 
that a BPS one-particle state is generically a bound state consisting of more than one 
charge centers, which are spatially distributed according to balance of classical forces 
[6]. Although initial studies were done for 1/4 BPS states in A^ = 4 theories and in 
classical setting, this was quickly elevated to a quantum statement [7, 8] and then 
extended to A^ = 2 theories [9, 10, 11]. From this new viewpoint, the wall-crossing 
occurs simply because the size of such bound states become infinitely large as a 
marginal stability wall is approached. Across the wall, these charge centers become 
mutually repulsive, precluding any bound state [7, 9].*^ Thus, what was initially 
thought to be a problem of " decay" is actually a problem of bound state formation, 
if viewed backward. 

For the simplest SU{2) example, the infinite tower of dyons and the massive 
vector meson which exists in the weak coupling side only should be all realized as 
loose bound states of the two BPS particles inside, monopoles and dyons of unit 
electric charge. At least this must be the right picture just outside the wall. 

The most interesting aspect of this problem is that rules for the formation or the 
dissociation of the bound state, whichever way one view it, seems quite universal and 
does not depend on detailed dynamics, which became more clear when the multi- 
center aspect was rediscovered later in the N = 2 black hole context [13, 14]. With 
N = 2 supersymmetry, BPS states are characterized by the charge 7, the central 
charge Z^, and the supermultiplet. Thanks to the partially preserved supersymmetry, 
the multiplet structure has a reduced form, which can be written as 

[j] ® ([1/2] © [0] © [0]) (1.1) 



'^"'^ These development were, however, restricted to weakly coupled theories, even though it gener- 
alized greatly the previous decades of monopole/dyon studies in the fully supersymmetric setting 
of iV = 2 and iV = 4 Yang-Mills theories. Much of these findings, and their relation to the more 
conventional monopole dynamics from 1980's and 1990's, was summarized in a review article [12]. 



The index that counts degeneracy of such BPS state is the 2nd hehcity trace 

fi = -Itr {{-ir^{2J,r) = {-ir^{2j + 1) . (1.2) 

For example, when a state 71 + 72 disappears across a marginal stability wall, and 
dissociates into 71 and 72 on the other side, the indices of these three kind of BPS 
particles are known to obey a universal formula [15], 

^-(71 + 72) = (-l)l<^^'^^>l-^|(7i,72)|fi+(7i)fi+(72) , (1.3) 

where ± denote the two sides of the wall. The wall itself is naturally defined by the 
condition of vanishing binding energy that \Zi\ + IZ2I = |Zi + Z2I. 

This simplest wall-crossing formula has been studied in many examples, general- 
ized to the so-called semi-primitive cases for 71 -|- A;72 states [15], and most recently 
embedded into an algebraic reformulation by Konsevitch and Soibelmann [16], which 
in turn was explained in more physical basis [17, 18, 19, 20]. Despite successes of these, 
later and more comprehensive, developments, much of the literature remain math- 
ematical and sometimes, from physics viewpoint, opaque. More direct approaches 
to the problem, based on physical mechanism of the bound states and the dissocia- 
tion thereof, do exist but until very recently were confined to special examples and 
situation. 

In Ref. [21], a new approach to the low energy dynamics of dyons in generic 
N = 2 Seiberg-Witten theory was proposed. Assuming that bound states of interest 
are large, which is always true whenever the theory is near a wall of marginal stability, 
the authors showed how a A/" = 4 supersymmetric dynamics can be explicitly written 
from the special Kaehler data of the vacuum moduli space only. When applied in the 
limit of a single dynamical probe dyon in the presence of another (very massive) BPS 
state, the bound states can be constructed explicitly and counted, again confirming 
the above (so-called primitive) wall-crossing formula. It is abundantly clear that 
this method can be used for an arbitrary number and varieties of dyons, as well, 
as long as the proximity to a marginal stability wall is satisfied. In this note, we 
wish to set up dynamics of arbitrary number of dyons near such a wall, with A/" = 4 
supersymmetry,*^ and generate wall-crossing formula via index theorem. 

Among the universal formulations that were previously attempted, Denef 's quiver 
dynamics picture [14] gives a very similar picture in the so-called Coulomb phase 
description. A recent work by Manschot et.al. [22, 23]*^ fully took advantage of the 
latter pictures and wrote down a wall-crossing formula. At the end of the day, our 



"^^To avoid confusion, we will use A/" to count quantum mechanical supersymmetries and N to 
count those of D = 4 field theory. 
'^'^See Ref. [24] for a short review. 



computation will lead to the same final wall-crossing formula. Since we start from 
Seiberg-Witten theory and derive the wall-crossing formula from scratch, we offer 
several improvements. 

The first improvement concerns the question of what is the relevant index theorem. 
In the Denef's Coulomb phase approach, the most comprehensive studies to date 
involve a truncation of dynamics where one ends up with a geometric quantization 
problem on the classical moduli space of charge centers, which are typically compact. 
In this paper, we denote such moduli spaces for n centers as M.n- For two-center 
case, this manifold is always S"^. The Lagrangian has no kinetic term, but a minimal 
coupling to certain magnetic field induces a symplectic structure on the moduli space, 
making it a phase space. In turns out, however, the naive low energy dynamics on this 
classical moduli space on M.n end up with too many fermionic degrees of freedom. 
The anticipated and empirically correct answer, which is a Dirac index [25], results 
only if one can somehow remove half of the fermions. This deficiency has remained 
unresolved until now. 

In this note, we will explain why the naive truncation to A^„ was ill-motivated. 
It turns out that there is no separation of scales, and all 3n bosons and 4n fermions 
are of equal massgap. Instead, one can choose to reduce the index problem to A^„ 
by deforming the theory with supersymmetry partially broken. As long as there is 
one supersymmetry left unbroken and since the quantum mechanics has a gap, the 
index is left invariant under the deformation. At the end of the day, we will thus 
have provided an ab initio derivation of the anticipated Dirac index on A^„, for the 
first time. 

The second concerns the physical interpretation of certain rational invariants, 
defined and extensively used by Manschot et.al. [22], of the form 

where the sum is over divisors of 7. The expression naturally appears in other for- 
mulation of the wall-crossing, most notably in Konsevitch-Soibelmann. In the course 
of enumerating the bound states of bosonic or fermionic statistics, we will encounter 
Q{P)/p'^ as a universal effective degeneracy of p identical particles of charge /3. It 
appears as the multiplicative factor from the normal bundle as one computes contri- 
butions from a submanifold fixed by the permutation group of order p.*^ 

Along the way, our work also clarifies relation between the field theory indices, 
namely the second helicity trace and the protected spin character, and the quantum 



'^^This same numerical factor 1/p^ had appeared before in the context of the D-brane bound state 
problems of 1990's [26, 27], where identical nature of the D-branes were also of some importance. 



mechanical ones. Quantum mechanical index usually suffers from ambiguity over 
the definition of (—1)^. Usual index formulae relies on certain (mathematically) 
canonical choice of (—1)^. Retaining three bosonic coordinates per dyons allow us to 
inherit both the spatial rotation group, denoted by SU(2)l, and the R-symmetry of 
N = 2 field theory, SU{2)r. The supersymmetries belong to (2, 2) representation, so 
both (— 1)^'^3 Qf SU{2)l and (— 1)^^^ of SU(2)ji are chirality operators. The second 
helicity trace is then computed unambiguously by Tr(— l)^'^^. We in turn relate the 
latter to Tr(— 1)^^^ which turns out to be equivalent to the canonical choice leading 
to the usual Dirac index formula. This derives, for the first time, the well-known sign 
pre-factors in the wall-crossing formulae universally. In addition, we also explain why 
the protected spin character of the field theory is actually computed by equivariant 
index, by showing that the quantum mechanical "angular momentum" operator that 
appears in the latter is actually a diagonal sum, J3+/3, from the spacetime viewpoint. 

The paper is organized as follows. Section 2 reviews Ref. [21] and generalize the 
low energy dynamics to the case of arbitrary number of dynamical charge centers, 
and note the universal nature of the potential terms. Section 3 defines the index as a 
method of BPS bound state counting, and in particular makes contact with the field 
theory indices, commonly known as the 2nd helicity trace and its generalization known 
as the protected spin character. It turns out that the quantum mechanics found have 
SU{2)l X SU{2)fi R-symmetry, each of which defines chirality operators (—1)^"^^ and 
(— 1)^^3_ TJ-^g ggy theory index corresponds to the former, while mathematical index 
formulae are more directly related to the latter. We discuss a universal relationship 
between the two, and conjecture that all BPS bound states in our quantum mechanics 
are all SU{2)ji singlets. 

Section 4 sets up index theorem for this dynamics and show how reduction to the 
classical moduli manifold may be achieved. Here we show why the naive derivative 
truncation leading to the geometric quantization is unjustified by demonstrating that 
there is no natural separation of scales between classically massive directions and 
classically massless directions. The main point is that A^„ is of finite size, and the 
quantum gaps due to this are always equal to those along the classically massive 
directions. We show, nevertheless, how one can deform the theory while preserving 
the index, such that classically massive modes are decoupled from the evaluation 
of the index, at the cost of partially broken supersymmetry. We also observe that 
the reduction process keeps a diagonal subgroup SU{2)j, and identify the generator 
<J3 = J3 + h as the operator usually used for equivariant index computations. This 
way, we show that the equivariant index of quantum mechanics on A^„ actually 
computes the protected spin character of A^ = 2 field theory. 

After the derivation of Dirac index in section 4, we go on to evaluate in section 5 
the wall-crossing formula by taking into account the bosonic or the fermionic statis- 



tics. Projection operators are introduced for the purpose, and the index formula is 
decomposed into additive contributions from various fixed submanifolds associated 
with coincident identical particles. The reduced index problems on the fixed sub- 
manifolds appears in the full index with a universal degeneracy factor ~ l/p^, which 
arises from orbifolding action of the p-th order permutation group S{p). Summing up 
all relevant contributions, we find an expression identical to Manschot et.al.'s wall- 
crossing formula. Section 6 further supports computation in section 5, by studying 
general orbifold index theorem. We close with summary and comments in section 7. 



2 A/" = 4 Moduli Mechanics for n BPS Objects 

In Ref. [21], a general framework for deriving moduli dynamics of dyons of Seiberg- 
Witten theory was given under the assumption that one works in the field theory 
vacuum where the central charge are almost aligned in terms of the phases of the 
respective central charges; in other words, very near the marginal stability wall. This 
program was then carried out explicitly when one can treat only one dyon as dynam- 
ical, with other dyons as external objects. In this note, we wish to generalize this to 
arbitrary number of charge centers, be they field theory dyons or charged black holes. 
For this, all dyons should be treated as dynamical, and we will denote their charges as 
7a's. For the above derivation of one dynamical center, the proximity to a marginal 
stability wall played an essential role, allowing the nonrelativistic approximation and 
thus the moduli space approximation possible, so we need to retain this assumption. 

While the moduli dynamics should have A/" = 4 supersymmetry, as demanded by 
the BPS nature of the dyons, simple off-shell A/" = 4 descriptions fail to accommodate 
key interaction terms. Furthermore, as we will see in section 4 where we compute the 
supersymmetric index, it is more convenient to take one of the four supersymmetries, 
say Qi, and give up others. For these reasons, we employ the M = 1 superspace [28] 
where this supersymmetry, Q4, is manifest. We package 3n bosonic coordinates, x^", 
and 4n fermionic superpartners, ^^"' and A"^, as 

^Aa ^ ^Aa _ -Q^Aa^ ^A ^ -^A ^ -q^A ^ ^g.l) 

with n auxiliary field 6^'s. The supertranslation generator and the supercovariant 
derivatives are then, 

Q = de + lOdt, D = de- lOdt . (2.2) 



2.1 Two Centers 

The general structure of two dyon dynamics can be inferred from the results in 
Ref. [21]. The latter actually derived the effective action of a single dynamical dyon in 
the background of an infinitely heavy core BPS state. When the core state consists of 
a single dyon, the effective action derived there can also be regarded as the interacting 
"relative" part C^^ of a two-dyon effective action, upon the usual decomposition, 

C = £"™- + C'^ , (2.3) 

where the trivial center of mass part was understood to be 

with Mtotai — ^ c>o understood. Here, let us recall basic structures of C^^ as dictated 
by the supersymmetry. 

£rei involves only three bosonic coordinates and four fermionic ones and can be 
further decomposed as 

^rel ^ ^el ^ £re« ^ (2.5) 

where 

with a = 1, 2, 3, and 

C{i = I de (i/C($)A - iW($)a/^$") , (2.7) 

with the condition 

dale = eab^dbWc (2.8) 

imposed. Note that this also implies dadalC = 0, which is solved by 

/C = /C(oo)-^. (2.9) 

\x\ 

We will see shortly how /C(oo) and q can be read off from the underlying Seiberg- 
Witten theory. 

As was claimed, this Lagrangian is invariant under four supersymmetries, 

5,b = -le^i^^ , (2.10) 



with four Grassman parameters e"^ and with ip'^ = A. The J\f = 1 superspace we 
employed is related to €4, so Cq and £1 are manifestly and individually invariant under 
these supersymmetry transformation rules. A less obvious fact, which is nevertheless 
true, is that the two are also individually invariant under all four 



dtCo^^ =0 = 6, 



I dtC{\ (2.11) 



if the auxiliary field h is kept off-shell. This is the feature that allows an easy gener- 
alization to n dynamical centers. The auxiliary field h takes the on-shell value, 

"- '^onshcll = T f^ + -XAfi'^'i''] ' (2-12) 



/ V 4 



which generates bosonic potential terms of type K? /2f and mixes up terms in £o,i- 
Nevertheless, J\f = 4 supersymmetry of £ = £0 + -^i still holds, now in far more 
complicated on-shell form. 



2.2 Seiberg-Witten 

Before we extend this to n dynamical dyons, we need to understand the role of the 
core-probe approximation and how it computes /, /C and W [21] in terms of the 
quantities that appear in the Seiberg-Witten theory. 

Let us consider a collection of charges '-/a, and represent it as a semiclassical 
state. The basic information about the semiclassical dyon state comes from the BPS 
equations of the Seiberg-Witten theory [29, 30, 31, 32] 

F,-<-iV0. = O, F;,-<-iWd = 0, (2.13) 

where F = B + iE with magnetic field -B's and electric field E^s, 0's are unbroken 
part of the complex adjoint scalars, each of which are labeled by the Cartan index 
i = 1,2, . . . ,r. -Fd's are defined through the low energy f/(l) coupling matrix as 

n - r^'F, , r^^ = -^ . (2.14) 

The pure phase factor (^ is determined by the supersymmetry left unbroken by the 
charge 7 in a given vacuum, and equals the phase factor of the central charge Z^ of 
the configuration. 

In a core-probe approximation, we split 7t = 7/i + X^a' 7-4' ^"^^ treat the latter 
n — 1 as a fixed background of total charge 7c = X^^/ 7yi'- As we saw in the previous 



section, the Lagrangian for the dynamical dyon (of charge 7/1) is characterized by 
three objects. 

The first is the mass function / = | Z^^^ \ as in 



1 „ ( dx^ 



(•- 



C--,f[^)^-, (2.15) 

where 

2,.=lt-<f> + lh-<f>D, (2.16) 

with the electric part 7^ and the magnetic part 7™ of the charge vector 7^^. The 
scalar fields here solve the above BPS equation with the other n — 1 charges 7yi''s as 
the background point-like sources. The fact we treat such dyons as point-like objects 
is justified by going very near a marginal stability wall, since this tends to separate 
charge centers far apart from one another. As we will see shortly, this proximity to 
marginal stability wall plays a central role in allowing us to construct nonrelativistic 
low energy dynamics of dyons. 

Clearly \ZyJ acts as the inertia of the probe dyon, which is position-dependent 
because of the background: this sort of identification is in accordance with general 
spirit of how one describe well-separated charged objects [33], which has been tested 
and used successfully for many soliton systems and even lead to exact moduli space 
metric in some cases [34, 12]. We also use the notation Z^ for the central charge of 
the charge 7 so that Z^ = Z^{oo). 

The other two, more important for the discussion of BPS bound states, are the 
potential K? /2f and the vector potential W, so that 

1 fdx\'^ K? dx - 

where these two are determined entirely by the charge distribution of 7a''s as [21] 

dW = *dlC, /C = Im[C-%J = Im[r^Zj-^— ^^, (2.18) 

with#5 

QhA' = (7h,7A')/2 



'^^This convention for the Schwinger product here follows the one used by Denef in Ref. [14, 15]. 
The original derivation of dyon dynamics from Seiberg-Witten theory in Ref. [21] used a different 
convention, such that 

(7,7') = (7,7')Dcncf = 2(7',7)Lcc-Yi- 

The tilde emphasizes the fact that the latter also used half-integral electric charges as opposed to 
integral ones, which is natural when we compute Coulomb energy. Magnetic charges are integral in 
either convention. 



for the Schwinger product. 

These are direct consequences of the equations (2.13), combined with the extra as- 
sumption of being near the marginal stabihty wall. Generically, the bosonic potential 
would have been 

\Zn\-MC^Zh] , (2.19) 

but this reduces to 

ie/2\Zh\ = {lm^[C^Zn]f/2\Zh\ , (2.20) 

as we move near the marginal stability wall defined by alignment of Z^ and Z^ [21]. 
The reason why we need this proximity to the marginal stability wall is clearly not 
because of inherent properties of the system, but rather because of the non-relativistic 
quantum mechanics approximation we employed. Far away from the wall, the po- 
tential energy would be not small compared to rest mass of the particles involved, 
which will bring dynamics to a relativistic one. However, we do not know how to 
handle interacting and relativistic particles at mechanical level. *^ Nevertheless, this 
approximation is good enough since we already know that BPS states are stable far 
away from marginal stability walls. 

An important subtlety we wish to point out here is the choice of C- In the core- 
probe limit, it appears that C = Cc = Z^^/\Z^^\ is the right choice, since we are treating 
7/i as an external particle in the background given by 7c = XIa'7^'- However, ^ is 
tied to the supersymmetry left unbroken by the configuration and furthermore we 
are interested in the supersymmetric bound states of 7c and 7/1. Around such a state, 
the low energy dynamics should have supersymmetries associated with 7t = 7c + Ih 
rather than those associated with 7c. 

One can understand this as capturing the backreaction of the background due to 
the probe. Failing to do so clearly will give us nonsensical answers since, otherwise, 
the supersymmetry of the bound state in question would not be aligned with the 
supersymmetry of the moduli dynamics. In the core-probe approximation, the two 
happen to be the same, Ct = ■^7t/I-^7tI — Cc simply because the total central charge 
is dominated by that of the infinitely heavy core state. As we give up the core-probe 
dichotomy, this accidental identity will no longer hold, and the preceding discussion 
tells us that one must always use Ct- 

As we give up the core-probe approximation and treat all charge centers on equal 
footing, the moduli dynamics will become quite complicated. The part of the above 
action that remains least affected by this extension is the Lorentz force, coming from 
—X ■ W type couplings. The coefficient g in W keeps track of how one particle's 
quantized electric (magnetic) charges see the other particle's quantized magnetic 



'^^Importance of the wall in the derivation of low energy dynamics of dyons was also recognized 
by others [35]. 
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(electric one) charges. W is Dirac-quantized and topological, and furthermore can 
arise only from sum of two-body interactions. Therefore, this part of the interaction 
can be reliably computed by adding up all pair-wise Lorentz forces, giving us 

with 

W^a = Yl ^^fiWf ™'=(f A - xb) , (2.22) 

where qab = {ia, 1b) f^ and yy^'^™'^ is the Wu-Yang vector potential [36] of a in flux 
Dirac monopole. Note that the Att flux of yy^*''"'^ dovetails nicely with half-integer- 
quantized qab, as demanded by the Dirac quantization. 

For general n also. A/" = 4 supersymmetry constrains the Lagrangian greatly and, 
as we will see shortly, the potential energy is tied to such minimal couplings. Knowing 
the latter will allow us to fix, almost completely, the analog of K? /2f as well. We will 
presently see how this works in n center case. A more difficult question is how the 
kinetic terms would generalize, to which we will only give a general statement rather 
than precise solution. In this note, our primary interest is in the supersymmetric 
index for non-threshold bound states, which is independent of details of kinetic term. 

2.3 Many Centers 

For many centers, it is more convenient not to separate out the center of mass coor- 
dinate. Let us label the centers by A = 1, 2, . . . ,n and denote their B? position as 
x^"" and the charge 7a. The A/" = 1 superfield content is 

^Aa ^ ^Aa _ -g^Aa ^ ^A ^ -^A ^ -q^A ^ ^2.23) 

with A = 1,2, . . . ,n and a = 1, 2, 3. A/" = 4 transformation rules are. 



6M = -le^V^ , (2.24) 

where as before '^^'^ = A"^. We again split the Lagrangian into the kinetic part and 
the potential part, 

£ = £o + /^i , (2.25) 

and look for £o i separately, with off-shell ft'^'s. 
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The n-center version of £i is, given (2.22), quite obvious, 

C, = fd9 (2/C^(<I>)A^-2Wa,($)D$^'^) , (2.26) 

since the second term gives precisely the Lorentz force among dyons and while the 
first is induced from the second by A/" = 4 supersymmetry; One can check easily that 

6, dtCi = (2.27) 

under all four supersymmtries, provided that 

dAalCB = leabcidAby^Bc-dBcy^At) , (2.28) 



and 



^abcdAbdscK^C — , 

dAadsalCc = 0, (2.29) 



for any A, B, C. We already learned that 



WAa = J2^^^^a'''''(^^-^'') 



B 

SO /C's also follow immediately via the A/" = 4 constraints as 

(7a, 7b) 



JCa = JCaM - i Y: M^ ■ (2-30) 

B ' ' 

Note that this obeys the constraints except at the submanifold, say A = {x^"- : 
xa = xb, {'Ja, 1b) 7^ 0}. The quantum mechanics can be very singular at such places 
also, meaning that we should excise A from i?^" and impose the regular boundary 
condition instead. 

It remains for us to determine /C^(oo)'s. These /C's and W's can be traced back 
to the original BPS equations (2.13), and found by keeping track of how motion of 
each center is affected by the presence of the other n — 1 centers. After solving the 
BPS equations, similarly as in the core-probe limit, we learn that 

ICa = Im [C'Za] = Im [C'Za] - ^ Yl r?^^ > (2-31) 

where Za is computed from the solution to (2.13) with the other n—1 charge centers 
taken as the background but, nevertheless, with the phase of the total charge, ( = 
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J2a ^a/\ J2a ^^I' ^^^^ ill the equations. As we noted above, this is because we must 
make sure to use the supersymmetries that are preserved by the bound state of all 
centers. This can be also seen from )Ca{oo) = lm.[(^^ZA], which allows J2a ^^i^ = 
as demanded by the antisymmetric Schwinger product. Note that this consistency 
condition would have been violated if we had used different C's for different /Ca's- 

The other piece Cq, containing kinetic terms, is a little more involved. The sim- 
plest way to find the most general Cq is via an A/" = 4 superspace. For this, note 
that the collection {$", A} can be thought of as dimensional reduction of a D = 4 
A^ = 1 vector superfield [37, 38].*^ In this map, x"'s come from the spatial part of 
the vector field, the fermions from the gaugino, and the auxiliary field b from that of 
the D = 4 N = 1 vector superfield. See appendix A for more detail. Here, we are 
mainly interested in A/" = 1 form of such a general Cq, which is available in Maloney 
et.al. [28], 

£o = J d9'-gAaBbD<^^''dt^'"'-lhABA^DA''-tkAaB^^''A'' + --- (2.32) 

where the ellipsis denotes four cubic terms that we omit here for the sake of simplicity. 
This Cq is also invariant under the four supersymmetries we listed above, 

6, I dtCo = (2.33) 

on its own with b^^s off-shell, provided that various coefficient functions derive from 
a single real function L{x) of 3n variables as 

QAaBti^) = [6:6l + e,\e'{)dAedBfm): 

hAB^ = S'^'dAadBbm) , 
kAaBi^) = e'^f^dAedBfm) , 

: (2.34) 

The A/" = 4 supersymmetry requires all those terms as well. See Appendix A for the 
complete form of Cq. 

Figuring out the precise form of L for n charge centers requires further work. 
For a single dynamical dyon in the core-probe limit, we know that it is related to 
the central charge function as d^L = \Z\. We expect that there exists a similarly 
intuitive generalization for n particles case as well. In this note, we are primarily 



"^^In this version of A/" = 4 superspace, L\ is not obvious. On the other hand, a more extended 
harmonic superspace form has been found to accommodate both kinetic terms and potential terms 

[39]. 
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interested in counting nonthreshold bound state, for which details of L does not 
enter. Determination of L can become an important issue, when we begin to consider 
non-primitive charge states. See next subsection for related comments. 

Again, the main point here is that Cq and £i are invariant under the four super- 
symmetries separately when we keep the auxiliary fields fc'^'s off-shell. Combining the 
two, it follows that the full Lagrangian 

Co + £i (2.35) 

is also invariant under all four supersymmetries. Integrating out fe'^'s generates po- 
tentials of type ~ K? and mixes up terms in Cq and £i, but A/" = 4 supersymmetries 
of the entire Lagrangian remain intact. 



2.4 Kinetic Function L : BPS Dyons vs BPS Black Holes 

Note that the potential part £i of the Lagrangian looks identical to the similar ex- 
pression previously found by Denef [14], which has been later used extensively for 
counting BPS black holes bound states [25, 22]. The latter relied on A/" = 4 quantum 
mechanical supersymmetry. Although we started with Seiberg-Witten theory for the 
derivation of £i, this part of Lagrangian is entirely determined by A/" = 4 super- 
symmetry combined with long-distance Lorentz forces among charge centers. Thus 
appearance of the same £i is hardly surprising. In fact, when we apply £i to BPS 
black holes, it is even more trustworthy, since the Abelian approximation that would 
underlie such an interaction form is valid all the way to horizon. One cannot say the 
same for field theory dyons, since at short distance non- Abelian nature must be taken 
into account. Nevertheless, as long as we are near a marginal stability wall and only 
long-distance physics matters, it is clear that £i is capable of describing both dyons 
and black holes. 

This does not mean that the moduli dynamics of BPS dyons and those of BPS 
black holes are identical. The difference resides in the kinetic part Cq of the La- 
grangian. As demanded by A/" = 4 supersymmetry, Cq is determined by a single 
scalar function L of the n position vectors xa- For instance, L for many BPS black 
holes of an identical charge was found by Maloney et. al. [28] 

L(fi, f2, • . . ) = -T^ I dx'^^ (2.36) 

iOTT J 

where ip = 1 + X^^(^/| xa — x\) with the mass m. On the other hand, for two- 
center dyon case, we expect smooth behavior near r*= [21] since, when the mutual 
distance is small, non- Abelian cores cannot be ignored and will smooth out Coulombic 
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singularities. Even if we use the naive Abelian results, d^L ~ 1/r at most. Comparing 
this to the two-body case of the supergravity result shows a substantial difference 
when the two objects begin to overlap. 

Indeed, there are situations when the two theories are expected to give different 
answers. There is no known example of A^ = 2 field theory dyon which is a bound 
state of two or more identical dyons. For black holes, however, no such restriction 
seems to exist. If a BPS black hole of charge 7 exist, we expect BPS black holes of 
charge A^7 also to exist, in fact with large entropy. In the present context of moduli 
quantum mechanics, the latter corresponds to a collection of many charge centers with 
many flat directions extending to spatial infinities and may be realized as threshold 
bound states thereof. In such cases, the kinetic term of the effective action at both 
short distances and long distances could be important. This problem is an important 
outstanding issue in wall-crossing phenomena in general, for it provided much-needed 
input data on what dyons or black holes are available, to begin with, to form bound 
states. 

Explicit forms of L for n BPS dyons and for n BPS black holes, respectively, will 
be studied in a separate work. 



3 R-Symmetry, Chirality Operators, and Indices 

We wish to compute index of the preceding quantum mechanics 

Tr ((-l)^e-^^) . (3.1) 

Since the quantum mechanics is gapped, of which much discussion will follow in 
next section, this quantity is truly independent of the parameter s. Thus, following 
the standard arguments, we will compute this in small s limit. Before proceeding, 
however, it is important to clarify what we mean by the operator (—1)''^. In order for 
the index to make sense, this operator needs to anticommute with supercharge(s), 

{(-ir,g} = o, 

which is the condition needed for 1-1 matching and thus cancelation between bosonic 
and fermionic states for nonzero energy eigenvalues. Clearly this is not enough to 
fix the overall sign of (—1)'^ on the Hilbert space, and an index is also plagued by 
this ambiguity. When we compute an index of standard Dirac operator or de Rham 
operators, there is usually a canonical choice that is used widely. We will come back 
to this, later in next section, but the choice is a matter of convenience only and, a 
priori, has no physical significance. 
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At field theory level, however, we have an unambiguous and useful definition of 
such an index, say, the second helicity trace, 

n = -^Tr{i-lf^-\2JsY) , (3.2) 

where the trace Tr is over a single particle sector of a given charge. We wish to fix 
the sign of the quantum mechanical index, in accordance with this. Irreducible BPS 
multiplets, tensor products of half- hyper- multiplet and a spin j multiplet, have the 
index 

fi([j]®([l/2]©2[0])) =(-l)2^-(2j + l), (3.3) 

so often we also write, 

n = Tr{{-lf^-') , (3.4) 

with the factored-out half-hypermultiplet understood. This naturally reduces to the 
low energy dynamics of dyons, which then must correspond to an index defined with 
a chirality operator that acts exactly like (—1)^"^^ 

n ^ Tr((-l)2-^«) , (3.5) 

but of course we need to ask here how such an operator is realized in the quantum 
mechanics. 

As can be inferred from discussions in Ref. [21], the quantum mechanics of previous 
section are equipped with 5*0(4) = SU{2)l x SU(2)ji R-symmetry. This is easiest to 
see in how the fermion bilinear couplings to rf/C and rfW combine to give, 

-^Cn^Aa/CB^^™^^" (3.6) 

in the component form, where, as before, ■?^^™=i'2.3 _ ^Aa=i,2,3^ ^aa ^ •\^a^ ^^^ ^ jg 
the 't Hooft self-dual symbol. The above form is precise when the metric is fiat, but 
appropriately modified preserving 5*0(4) symmetry when it is not. For each particle 
indexed by A, bosonic coordinates are in (3, 1) representations while the fermions are 
in (2, 2). Since spatial rotations rotate x"^ as 3- vectors, SU(2)l should be interpreted 
as the rotation group, while SU(2)ji must be descendant of SU(2)ji R-symmetry of 
the underlying Seiberg-Witten theory. The latter rotates only fermions and leaves 
the position coordinate intact.*^ 

In particular, the four supersymmetries are labeled by the 5*0(4) vector index, 
and thus are in (2,2) representations. Denoting generators of these two 5'f/(2)'s by 
J and J, respectively, we thus find 

{i-ir\Q} = = {(-1)2^3^ g|. (3.7) 



*^In the core-probe approximation of Ref. [21], only SU{2)ji were generically there, but this was 
an artefact of treating some of dyon centers as fixed background. 
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The quantum mechanics have two unambiguous and physically meaningful chirality 
operators that can be used for index computation. The desired (— l)^"'^ is one of 
them, therefore, we have an unambiguous way of computing the field theory index 
from the low energy quantum mechanics. 

On the other hand, there is an interesting and universal relationship between these 
pair of chiral operators in the quantum mechanics. Restricting our attention to the 
relative part of the low energy dynamics again, we have 

(-1)2"'3 = (_X)Ea<s(7a,7b>+«-1(^_]^)2/3 _ (^3 3) 

where n is the number of dyons in the low energy dynamics. This is easy to see by 
considering how the two SU{2) generators are constructed in the quantum mechanics. 
For SU{2)fi, which rotate only fermions, we have 

^•^ = E (-^ ''^^^ [^^'' ^^1 + i [^^"' ^l) ' (3.9) 

where the hat signifies the unit normalized fermion. The spatial rotation generators 

Ja = L, + Yl (-^ '-^^ [^^'' ^""1 - i [^^"' ^^]) ' (3.10) 

are similar but differ in two aspects: first, since SU(2)l rotates xah, the generators 
include the orbital angular momentum L; secondly the fermions rotate differently, 
as reflected in the sign of the last term. This latter difference generates a relative 
sign between the two chiral operators for each (2, 2) representation of fermions, thus 
explaining (— 1)""^. The other sign is equally simple, and come from well-known 
piece of charge-monopole physics, where the orbital angular momenta is schematically 
something like 

L ~ > ajA X TT^ + > — ^- 3.11 

■^^ -^^ 2 \xa — xb\ 

A A>B ' ^ ■"' 

with the covariantized momenta tva- The orbital angular momentum is constructed 
from tensor product of spin (7a,7_b)/2 representations times usual integral angu- 
lar momentum. Then regardless of which particular SU{2)l multiplet the state is, 
integrality vs half-integrality of the orbital angular momentum is unambiguously de- 
termined as 

(-1)2^3 = (_l)i:A>i5(7A,7S> _ (3 12) 

Note that this does not require L being symmetry operators. 

Thus, we have the second helicity trace of A^ = 2 dyons which can be computed 
via the low energy quantum mechanics as 

n = TT ((- 1)2-^3 e-^ff) = (_i)EA<s(7A,7B>+n-i X Tr {{-ly^-'e''^) . (3.13) 
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In the subsequent computation, with this relation in mind, we will eventually identify 
(— 1)^^3 as ^j-^e canonical chirality operator (—1)-'^. For this, there is another sign issue 
to settle, later when we begin to quote index formula from literature, since the latter 
come with a canonical choice of (—1)^, which may or may not equal to our choice, 
(— 1)^^3^ ]3^^ -^g postpone this to end of next section. 

Another reason why (—1)^'^^ is useful, even though we ultimately want (—1)^"^^, 
can be found in the observation [41] that all explicitly constructed field theory BPS 
states, to date, are in SU{2)ji singlets times the universal half-hypermultiplet (from 
the center of mass part in quantum mechanics viewpoint). If this is generally true, we 
can see that the index with (— 1)^^^ jg always positive and truly counts the degeneracy. 
An interesting question, therefore, is whether in the low energy quantum mechanics 
we derived all supersymmetric bound states are SU{2)r singlets. 

An interesting variant of the second helicity trace is the protected spin character 

[41], #9 

Tr ((-1)2-^3^2 J.+2/3) ^ (3 14) 

where again we took out the universal half-hypermultiplet from the trace for simplic- 
ity. This clearly reduces to, in quantum mechanics, 

Tr ((-l)2'^3^2j3+2/3) _ (3,15) 

Later we will also see how this quantity is naturally computed, after we reduce the 
index problem to the more familiar one that relies only the classical moduli space 
/C = 0, by the equivariant index that counts "angular momentum" representations. 
As we will see, this reduction process cannot carry the entire A/" = 4 supersymmetry, 
and, of 5*0(4) R-symmetry, only a diagonal SU{2) subgroup generated by J + / 
survives as global symmetry. The equivariant index on /C = space does not count 
representations under spatial rotations but under simultaneous rotation of spatial 
SU{2)l and N = 2 R-symmetry SU {2) r."^^^ See section 4.4. for more detail. 



4 Index Theorem for Distinguishable Centers 

Now we turn to the problem of counting ground states of the above quantum me- 
chanics, or equivalently counting BPS bound states of n dyons. Since the quantum 



^^We are indebted to Boris Pioline and Jan Manschot for bringing the question of the protected 
spin character to our attention. 

^^°0f course, if the SU(2)ji singlet hypothesis actually holds for the ground state sector, the end 
result would not know about /s, anyway. In fact, on the basis of this hypothesis, this equivalence 
was anticipated previously [22]. Our argument in section 4.4 will prove the identity without such 
an assumption. 
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mechanics has a potential, ~ /C^, one may expect that the problem can be reduced 
naturally to another problem on the classical moduli space of 2(n — 1) dimensions, 
say, 

Mn = {a:^" I /C^ = 0, A = 1, 2, . . . , n}/R^ , (4.1) 

where the division by R^ is to remove the flat center of mass part. This classical 
moduli space is generically a little more complicated since some of the centers could 
be associated with identical particles, which we will deal with in the next section. 

This reduction is not as straightforward as one might think, however. Ref. [25], 
for example, suggested that one can ignore the (then unknown) kinetic part of the 
Lagrangian. Effectively, in our notation, this would involve a geometric quantization 
of £i, 

which is obtained by truncating higher-derivative parts in Cq. The auxiliary fields, 
fe'^'s, are now Lagrange multipliers, imposing ICa = as constraints and leaving a 
lowest Landau level problem on A^„ with the magnetic fields XIa'^^-4- However, 
computation of the resulting index, if we take Cgeometric verbatim, generates wrong 
results relative to other known spectrum; The geometric quantization of Cgeometric 
would lead to index formula that is known to generate empirically incorrect answers. 

For two body case, for example, the degeneracy 2|g| has been known to be the 
correct answer for many explicit constructions. See, for example, Ref. [9] and Ref. [21] 
for explicit two-dyon bound state construction in the weakly coupled and in the 
strongly coupled regions of Seiberg-Witten theory, respectively. On the other hand, 
the naive lowest level Landau problem (or equivalently the geometric quantization 
problem) gives 2\q\ + 1. One would hope that the effect of fermions in Cgeometric will 
fix this, but this apparently does not happen. 

The truncation of the kinetic terms in the presence of fermions is quite subtle, 
since while bosons acquire a symplectic structure thanks to the magnetic field, there 
is no such analog for fermions. Setting the kinetic term of fermions to zero will cause 
the canonical commutator ill-defined, making the whole reduction process ambiguous. 
One can try to reinstate kinetic terms on A^„ as a regulator, but then, the main issue 
is that the number of fermions in Ci is 4(n — 1) real while the number of bosons is 
2(n — 1) real, and these lead to de Rham cohomology problem on M.n- For not too 
small q and when M.n is Kaehler, for example, the index of such a quantum mechanics 
coincides precisely with the state counting of the bosonic geometric quantization 
problem,*^^ again giving us wrong result for the index. 



"^^^Sce for example Ref. [40], where in effect a regularized version of these problems were considered 
with kinetic terms on A4„ and for its fermionic partners present. 
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Really at the heart of the problem is, however, the fact that the classical massive 
directions are in fact no more massive than the classically massless directions. Be- 
cause the classical moduli space A^„ is of finite size*^^, it comes with various gaps at 
quantum level, and it so happens that these quantum gaps are one-to-one matched 
and identical to the gaps associated with the classically massive directions: the dy- 
namics cannot be really split into two distinct sectors of heavy and light modes, at 
all, and contrary to initial expectation, the reduction to A^„ cannot be justified. 
In fact, this lack of separation of scales is easiest to see in how fermions enter the 
Hamiltonian. Half of fermions get mass from dK, while the other hand get mass from 
dW. However, A/" = 4 supersymmetry of the quantum mechanics tells us that the 
two are one and the same object, and fermions coupling to dK, are no more heavier 
than those coupling to dW. 

Fortunately, we can still decouple these classically massive directions in the com- 
putation of the index problem. This involves a deformation that breaks all but one 
supersymmetry, yet because the quantum mechanics is gapped and the surviving 
supercharge is effectively a Fredholm operator, it can be done while preserving the 
index. Later in the section and in Appendix B, we explicitly show that, as far as 
computation of the index goes, we may reduce the moduli dynamics to an effective 
A/" = 1 supersymmetric quantum mechanics with target Ain-, 

2 ^ 2 

where ^ is a gauge field on A/l„ such that 



^forindcx only (-Mn) = ^ G.^z'^z'' + ^ G,.rV + " " " " A,Z>^ + ■ ■ ■ , (4.3) 



dA = J^ = d\ > WAadx' 



(4.4) 



/ M„ 

and G is the induced metric on Ain- This Lagrangian must be used only for the 
purpose of computing index. 

The key point here is that the number of fermions is exactly half of that in 
j^geometric- Siuce these fermious live on the tangent bundle of A^„, we have a nonlinear 
sigma model with real fermions. The relevant wavefunctions are spinors on A^„ and 
the index in question becomes a Dirac index, 

M{ia}) = I Ch{J')A{Mn) = I Ch{r) (4.5) 

J Mn J Mn 

with the Chern character Ch of J-". A is the A-roof genus of the tangent bundle, 
which will be shown to be trivial for all A4„'s. This formula counts the index when 



^-'^^There are also some exotic cases corresponding to the scaling solutions. In these cases, the 
moduli space is non-compact, from short distance side, but its volume in the naive flat metric is 
still finite. 
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we view individual charge centers as distinguishable; in section 5, we will extend the 
formula appropriately when identical particles are involved and along the way see 
why the rational invariants of the form ~ Vt/p^ with integer p > 1 appears in various 
wall-crossing formulae. 

The Dirac index found here is consistent with de Boer et. al.'s observation [25] 
that empirically correct answers emerge for n = 2 and n = 3 if one assumes that the 
relevant quantum mechanics admit spinors on Ain as the wavefunction. This can be 
then generalized to the refined index (or equivariant index) and make contact with a 
series of recent works by Manschot et.al [22, 23]. 

4.1 Two Centers: Reduction to S'^ 

Supersymmetric ground states were found and counted for n = 2 case in Ref. [21], 
which gave the correct answer of 2g at the end of the day. As expected, the wave- 
functions are all maximized near the classical "true" moduli space /C = 0, which was 
nothing but a two-sphere threaded by a flux of Anq. However, the wavefunctions can 
also be seen to be very diffuse, too much so to let us call it "localized" there. 

Here, we will illustrate why a naive reduction to A^2 = "S*^ by throwing away 
entire kinetic term is wrong. After the latter procedure, one ends up with Cgeometric 
for which we need to either geometrically quantize over S*^ or regularize the dynamics 
by reinserting kinetic term on S*^ and concentrate on the lowest Landau level. If we 
follow the second viewpoint, we end up with a two dimensional nonlinear sigma model 
with four real fermions, so effectively we will have thrown away only the bosonic radial 
coordinate from the original moduli dynamics. 

Let us consider the zero point energy of the relative part of the two-center me- 
chanics. Three bosons can be split into "radial" directions, on which /C and the mass 
function / depend, and flat "angular" directions. With K = a — q/r and positive 
a and g, the ground state is at r = ro = g/a, and the radial direction becomes a 
harmonic oscillator of frequency w = a^//(ro)g, so 



'■'radial 



^rad.al + U > -^^ . (4.6) 

2/ 2/(ro)g 

The angular part, although classical flat, also comes with a gap due to the flnite 
volume, and the energy quantization there goes as 

f 2 _ 2 2 

E , ~ - ^ > —^ (A 7) 

-^sphere — n j: f \ 2 — c^ £ / \ i V*- ' 

2/(ro)r^ 2f(ro)q 

since the angular momentum is bounded below, in the presence of the flux, by q. The 
four real fermions are paired up into two fermionic oscillators of the same frequency 
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w as above, so we get contribution from the fermion sector as 

E fermion ^ ("^/ +m'f~l)w> - , (4.8) 

^ -^ ^ /(ro)g 

where we again see that there is only one scale in the fermion sector also. Of course, 
the behavior of fermionic degrees of freedom must be the same as the bosonic ones, 
since we have supersymmetry. 

This shows that, without further deformation, the gap of the classically massive 
radial direction is exactly the same as the rest of the degrees of freedom. If we wish to 
localize the problem to A^2 = S'^ by removing the radial mode, we must do something 
else so that the gap along the radial direction and the gap along A^2 are different, but 
this seems impossible under the A/" = 4 supersymmetry of the quantum mechanics. 

Let us remember here that, for the evaluation of index, one needs only two things: 
a Dirac operator of some kind and a chirality operator that anticommutes with it. 
One would like to compute the index 

Tr(-l)^e-*^ , (4.9) 

for interacting part of the theory. Let us, for the sake of definiteness, take H = Q^, 

and evaluate 

Tr(-l)^e-"'3l _ (4 10) 

Af = 4 supersymmetry is useful since it constrains dynamics but all of them are 
not really necessary to define an index. It is clear that, as long as we preserve this 
quantity, we can even break A/" = 4 supersymmetry. 

Of course. A/" = 4 supersymmetry is important when it comes to generating correct 
supermultiplet structure to the bound state, but that only concerns the free center 
of mass part. The index must be computed from relative interacting part of the 
dynamics, only for which we will break A/" = 4 supersymmetry. 

Thus we are motivated to give up rf/C = *dyV condition, thereby keeping only 
Q = Q4: unbroken. Let us replace 

/C^^/C (4.11) 

with some arbitrarily large number ^ while keeping W as it is. The ground state 
energy counting is now 

^w w C,w + w 

-^radial ~r ^sphere ~r -t-' fermion ■— „ "o" O ' V / 

since the half of the fermions (A and ■?/'^) get the mass from d{^JC) and the other half 
from dW. The angular momentum sector mass-gap, w/2 = q/2f{ro)rl, is unchanged 
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since the classical vacuum, /C = and thus the radius Tq, and W are intact under 
this rescaling. 

It is not difficult to see that the reduced dynamics, after integrating out heavy 
modes, is a A/" = 1 nonlinear sigma model onto M.2 = S'^ coupled to an external 
vector field W. See Appendix B for complete detail of the reduction process. We 
note that since Ai2 = S"^ happens to be Kaehler, the unbroken supersymmetry gets 
accidentally extended to A/" = 2, although this is not important for our purpose. 



4.2 Many Centers: Reduction to A4n 

Similarly, we wish to deform the theory by rescaling ICa — > ■C^a, when we have many 
dynamical charge centers, as well 

+^iK,Am^'' - ^W^a($)D$^"^ , (4.13) 



where ^ is an arbitrarily large number. As in the two-center case, the bosonic poten- 
tials are quadratic in /C^'s and there are n — 1 "radial" directions that are of mass 
~ ^. There are also 2(n — 1) fermions that couple to (i(^/CA)'s, so they are also of 
mass ~ ^. The two sets can be decoupled together, thereby reducing the index prob- 
lem io M-n with real fermions. It leaves behind a A/" = 1 supersymmetric quantum 
mechanics onto Mn with 2(n — 1) bosons and 2(n — 1) real fermions. The process 
does not affect the free center of mass part, so the latter still comes with 3 bosonic 
coordinates and 4 fermionic ones. 

We may further deform the kinetic part, £o, by taking the simplest form of the 
kinetic function, 

L = - ^ m^x^ ■ x^ , (4.14) 

^ A 

which amounts to 

9AaBb = SAB^abmA, flAB = ^ABmA, kAaB = 0, (4-15) 

and setting cubic terms to zero as well. The simplest way to justify this deformation is 
that the kinetic function approaches this flat metric when distances between charge 
centers approach inflnity. This asymptotic form is more than good enough since 
we can always tune the fleld theory vacuum, so that we stay arbitrarily near the 
marginal stability wall. There, Im[(^~^Z^] approaches zero, and the submanifold Ain 
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is arbitrarily large. Since the index cannot change under the continuous and sign- 
preserving deformation of Im[(^^^Z^], and since the ambient metric is effectively flat 
for large Ain, the index will be unaffected by this choice of metric. 

This leaves us with a very simple A/" = 1 quantum mechanics 

^deformed = J d9 (^mAD^^''dt<^^'' - ^TUaA'^DA^ 

+i^lCAmA''-iWAamD<^^A , (4.16) 

with target R^"" modulo submanifolds given by ICa = ±00. Of this, the free center 
of mass positions R^ and the accompanying four real fermions decouples, leaving 
behind the interacting part of the moduli dynamics onto R^^"'~^\ This free part is 
also essential since it generates the basic BPS multiplet structure (whose content 
equals half of a hypermultiplet) to the bound state. Then, by taking ,^ — )■ 00, we 
decouple n — 1 "radial" directions and 2(n — 1) accompanying heavy fermions, and 
end up with a nonlinear sigma model onto A^„ with real 2{n — 1) fermionic partners. 
See appendix for detailed derivation of this fact. 

Thus, we arrive at the effective Lagrangian, which can be used for the purpose of 
computing the index of the original n center problem, 

<rindexo„iy = ^ G,.z^z^ - A,z^ + '- G.^rv + \ G^.^i^F^^^^ + \ T.^rv 

(4.17) 
again with the induced metric G on A^„ and, as we already noted, 

(4.18) 



dA = J= = d\Y^ yvAadx^" ] 



A4„ 



Since each Wa is a sum of Dirac monopoles at if = xb's, we find 
^ = rf('$^$^^^^^^Wf™^(f^-fB)rfx^" 

V A B^A / Mr. 

= d[Y^ hAlM Wf ™^(f A - xb) dix^'^ - x^'^) 
\a>b 2 

= $^^^^^^^''^™^(:?a-:?b), (4.19) 

A>B 

where J^^^'^"-'^ is the Dirac monopole of flux An. Of four supercharges, (^4 survives the 
deformation process above, which is then further reduced to Qmu as heavy modes 
are integrated out. 
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4.3 Index for n Distinguishable Centers 

Since this is the plain old nonlinear sigma model twisted by the minimal coupling to 
A, the reduced supercharge is represented geometrically as the Dirac operator with 
a [/(I) gauge field 

Q4 ^ Qa^„=7'^(^V^ + A) , (4-20) 

whose index, according to Atiyah-Singer index theorem, is given by 

X„({7a}) = Tr ((-l)^-.e-^«') = / Ch{J^)A{Mn) , 

^ ^ JMn 

as promised, where we must assume a canonical choice of the chirality operator. This 



IS, 



[-ifM^ =(2i)"-l^l---^2(n-l) ^ 



in terms of properly normalized and ordered fermions. See next subsection for how 
this choice squares off with physically motivated chirality operators (— l)^"^^ and 
(— 1)^^3 of section 3 and how the latter chirality operators reduce on A^„. 

Curiously enough, the A-roof genus A does not contribute to the index, thanks 
to the simple topology of A^„. To see this, let us first note that the ambient space, 
in which A^„ is embedded is essentially R^"-. For instance, take Xi = to remove the 
translation invariance and make the ambient space i?^("~^\ and then impose ICa = 0, 
of which n — 1 are linearly independent. Therefore, Ain is a complete intersection in 
^3(n-i)_ gjj^cg A-roof genus is a multiplicative class, we have an identity, 

A{TMn)A{NMn) = A (Ti?='("-i)|_^J = 1 (4.21) 

among the tangent and the normal bundles. However, rf/C^'s are nowhere vanish- 
ing normal vectors on A^„, and thus the normal bundle NAin is also topologically 
trivial*""^^, and 

AiTMn) = 1 . (4.22) 

It is important to note that this decoupling depends only on the topology of the 
ambient space, namely the original 3n dimensional moduli space, near the surface 
Ka = 0.'^^'^ This triviality is also implicit in the explicit evaluation of these Dirac 



#i3We are indebted to Bumsig Kini for pointing this out to us. 

"^^"^Note that similar argument wih not lead to triviality of other multiplicative class since typically 
they require complex bundles in order to be defined. For instance, Td{A42k) or c(A^2fe) cannot be 
argued to be trivial in this manner, for instance, since the normal bundle of A42k inside the relative 
position space E?y^^~^> is of odd dimension and, if irreducible, cannot be complex. In particular 
M.2 = S^ has a real line as the fibre when embedded into B? , which is consistent with the nontrivial 

Cl- 
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indices in Ref. [23]. Although it turned out there that A(TA^„) factor did make an 
important difference for evaluation of the equivariant index, the non-equivariant limit 
is consistent with trivial A. 



4.4 Reduced Symmetry, Index, and Internal Degeneracy 

Since we arrived at the nonlinear sigma- model on A^„ only after the deformation 
of the dynamics, which in particular removes the extended supersymmetry, we must 
first ask whether various operators survive this procedure of deformation and the 
subsequent reduction process ^ — t- oo. Of the four original supersymmetries, Q4 
survives the deformation. It's on-shell form will be smoothly deformed as well, which 
goes like 

Q^ = ... + A^/C^ -^ Q^ = ... + ^A^/C^ . (4.23) 

The ellipsis denotes parts unaffected by the deformation. We emphasize again that 
this supersymmetry is explicitly preserved since the deformed Lagrangian (4.13) is 
written in the superspace associated with (^4. Then, given that Q4 is a gapped elliptic 
operator, at ,^ = 1, this deformation preserves the index as we increase C, [11]. This (^4 
reduces to Qm^ of the nonlinear sigma model on Ain, and obviously the Hamiltonian, 
Ql/2, gets similarly deformed and eventually reduced to the natural one on A^„. 

This leaves the global symmetry operators and the chirality operators. With the 
A/" = 4 supersymmetry partially broken, the 5*0(4) R-symmetry can be easily seen 
to be broken. On the other hand, the deformation commutes with rotation of Xyi's, 
so we expect to see some SU{2) symmetry does survive the process. The question 
is which SU{2) in 5*0(4) = SU{2)l x SU{2)ji remains unbroken. The answer is the 
diagonal subgroup, SU{2)j, generated by 

Ja = Ja + Ia. (4.24) 

One can see this in several different ways. 

Firstly, both J {SU{2)l) and / {SU{2)fi) are broken by themselves, since they 
both act nontrivially on heavy fermion sector. The diagonal generators jT^'s, on the 
other hand does not involve A fermions and leave the heavy sector ground state 
untouched. Secondly, after deformation and reduction to Ain, the dynamics is a 
nonlinear sigma model, where fermions transform identically to bosons. Recall that 
bosons and fermions used to belong to (3,1) and (2,2) of SU{2)l x SU{2)r. In 
the reduced dynamics, symmetry properties of the bosons and fermions cannot be 
different, and indeed under the diagonal subgroup, bosons and fermions transform 
identically. Finally, after the deformation, the dynamics has only one real supersym- 
metry (^4 so no R-symmetry is expected. However this supercharge originates from 
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a (2,2) multiplet under SU{2)l x SU{2)ji, so has to transform nontrivially under 
either of the two individually. On the other hand, becaues JT does not rotate A's, JT 
commutes with Q4 and also with its reduced version Qm„- At the level of reduced 
dynamics on Ain, this SU{2)j is not an R-symmetry but a global symmetry that 
arises from the universal isometry of Ain- 

While we are on the question of symmetry, let us digress a little and consider the 
equivariant index or refined index one encounters in literature on wall- crossing, of the 
generic form 

Tr ((-1) V) (4.25) 

with a "rotation" generator js along z-axis. Most such computations are based on 
some version of low energy quantum mechanics on the classical moduli spaces, our 
MnS, but as we saw above, the "rotational symmetry" of A^„ is in fact not the purely 
spatial rotation but a diagonal subgroup of spatial rotation SU{2)l and the field 
theory R-symmetry SU{2)fi. Therefore, the refined indices that have been computed 
are in fact 

Tr ((-1) V^O = Tr ((-1) V'''+''') (4.26) 

so actually would equal the protected spin character 

Tr((-l)'-^3^2J3+2/3) (4 27) 

of the field theory, if we are allowed to choose the chirality operator (—1)'^ of the 
quantum mechanics to be (— l)^*^^. 

So this brings the question of what happens to the two natural chirality operators, 
(_X^2J3 ^^^ ^_iyh^ when we deform and reduce the dynamics in favor of a A^„ 
nonlinear sigma-model. As we saw, the two SU{2) symmetries are lost individually, 
so operators like J3 and I3 can no longer be used to classify eigenstates. Nevertheless, 
(_X)2'^3 g^Yid (—1)^'^^ are still sensible chirality operators. Even after the deformation, 
one can show directly (—1)^'^^ as a product of all fermions while (— 1)^"^^ is again the 
same product of all fermions times (— 1)Sa>s(7a,7s>+"--i^ Both anticommute with 
the surviving supercharge (^4, so still defines chirality operators. This is not much 
of surprise since they simply measure the most rudimentary information about the 
states, i.e., whether, before deformation, the state was in a integral or in a half-integral 
representations. 

When we reduced the dynamics to Ain, however, we must properly redefine these 
chirality operators by evaluating them on vacuum of the heavy oscillators. For in- 
stance, consider (—1)^'^^ for the simplest n = 2 case. The canonical chirality operator 
on JUn = S"^ is, as noted before. 
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with the natural orientation arising from embedding of S^ to R^. To relate this to 
(— 1)^^3^ we remember to set the heavy fermions, ip^ and A, to their ground state, 
which gives precisely 

it turns out.*^^ Clearly, we may repeat this for each sector of 4 fermions labeled by 
A, and find 

m-irmhea.y = l[^^^A^\ = ("1)^"^" • (4-28) 

A 

Therefore, the chirality operator (— l)^'^-^ prior to the deformation, smoothly descend 
to the canonical chirality operator on ( — 1)^-^", upon deformation and subsequent 
reduction of dynamics, and leads to the standard Dirac index X„. 

Since the desired index Q needs (—1)^"'^ as the chirality operator, we then use 
(3.13) to relate (— 1)^"^^ to ( — 1)^^^, and find an unambiguous answer, 

^distinct ^ (_l)EA>B(7A,7s>+n-l ^ X„({7a}) . (4.29) 

On the left hand side, we emphasized the fact we are yet to incorporate the statistics 
issue. We will see in next section how this generalizes when we impose statistics to 
the index computation. Before asking the question of statistics, however, there is 
still one more ambiguity to the expression above, since so far we did not take into 
account of the internal degeneracy and quantum numbers of the individual charge 
centers. The left hand side is still defined with respect to (—1)^"^^, so adding internal 
degeneracy factor can be accomodated by writing 

n 
^distinct ^ (_l)EA>s(7A,7s>+n-l x X„({7^}) X JJ fi^ (4.30) 

A=l 

where individual l^^'s are also computed as the trace of ( — 1)^"^^ (as usual modulo the 
universal half-hypermultiplet part). As usual, we assume that there is no significant 
coupling of these internal degeneracy to the quantum mechanical degrees of freedom. 
Sometimes, we will also write this as 

^_^distinct^ X (-l)^A>B(yAnB) = X„({7^}) X JJ(-fiA) (4.31) 

A 

which is more convenient when keeping track of statistics, since, for SU{2)ji singlets, 
the Bose/Fermi statistics are naturally correlated with the sign of — fi^'s. 



^^^This can be seen most easily when we choose the ordering of ja^s such that {jaj'^b) are ah 
nonnegative for A > B, which is also the convention chosen in Ref. [22] for non-scaling cases. 
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5 Index with Bose/Fermi Statistics and Rational 
Invariants Q 

So far, we pretended that dyons involved are all distinct, and studied the supersym- 
metric bound states thereof. In reality, this is not quite good enough since we often 
need to understand bound states of many identical dyons, obeying either fermionic or 
bosonic statistics. The effective true moduli space, for example, has to be an orbifold 
of type 

Mn/T (5.1) 

where F is a union of permutation groups that mix up labels for identical particles, 
with proper action on wavefunctions. Equivalently, the index should be computed 
with appropriate projection operator inserted, 

fi = Tr {{-lf^^'e"^Vr) (5.2) 

where Vr projects to wavefunctions obeying either Bose or Fermi statistics under the 
exchange of identical particles. 

The orbifolding reduces the volume of the moduli space, so given the index formula 
which is an integral over the manifold, we should expect to see factors like l/d\ as 
a result of having d identical centers. However, action of F is not everywhere free, 
since when identical particles are on top of one another, the action is trivial. There 
are complicated fixed submanifolds under F, making the problem very involved, and 
in particular there should be additional contributions from the fixed manifolds under 
the orbifolding action. 

5.1 The MPS Formula 

Before we carry out such a computation directly, it is instructive to recall a recent 
result by Manschot, Pioline, Sen (MPS), who evaded this complication algebraically, 
and replaced it by a sum of many index problems with distinct charge centers [22]. 
They argued that one can recover the correct index, by adding indices for a series of 
artificial problems with a smaller number of charge centers. In this set of effective 
index problems, the trick requires the following rules: When the reduced problem 
has d particles of the same kind, MPS divides the index by l/d\. When one has a 
particle of nonprimitive charge as a part of such a reduced problem, one must also 
use, in place of the true intrinsic degeneracies Q of the particles, a mathematical one 

n(-r)^j:^^^, (5.3) 

pn 
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where the sum is over the positive integer p such that 7/p belongs to the quantized 
charge lattice of the theory. Note that ^(7) = ^(7) whenever 7 is primitive and 
^{Pl) = ^{l)/p'^ if 110 non-primitive charge state exists. 

For illustration, let us take two primitive charges (3i and (32- Suppose that, among 
all possible linear combinations of the two, only these two states exist on one side of 
the marginal stability wall. Labeling the degeneracy by ± depending on which side 
of the wall we are considering, we thus assume that 

fi+(m/3i + A;/32) = 0, unless (m. A;) = ±(1, 0) or (m. A;) = ±(0, 1) . (5.4) 

The sign of r2i^2 = ^^(/3i,2) are correlated with the statistics assignment of the 
particle; a hypermultiplet has Q = 1 and must be treated as fermions while a vector 
multiplet has Q = —2 and must be treated as bosons. Under this assumption, we 
have Q{pPi^2) = ^{Pi,2)/p'^- Manschot et.al.'s formula then simplifies to, 

-Q-{mf3i + k(32) X (-1)Sa>b(7a,7s> 
= 4-X„+,(/3i, /3i, . . . , /32, /32, ■ ■ ■ ) (-f^i)™(-^2)' 

+ 7 Wi ^rn-l+k{2Pu /3l, /3l, . . . , /32, /32, . . . ) ^ {-n^r~'{-n2)' 

+ 7 \wrr; X„.-2+fc(3/3i, /3i, /3i, . . . , /32, /32, • • • ) ^ (-l]i)"-3(-fi2)' 

[m — 6)\K\ 6'^ 

+ ^y^^^^^^y|^X„.-2+fe(2/3i,2/3i,/3i,/3i...,/32,/32,...) (^) {-n,r-\-n2f 

+ ••• , (5.5) 

where the sum is over all unordered partitions of m(5i + k(52 respectively, although 
we listed above only part of the partitions of m. For the overall sign, we re-labeled 
the individual charges /9i, /9i, . . . , /92, /32, ■ ■ ■ and called them 7^'s. This sum and each 
term in it can be characterized by the following set of rules: 

(i) The sum is over all unordered partition of m/3i + A;/32 = Xls ^sPs where Ps = 
{PsiPi +Ps2P2)- For each P^, we will have a factor 0,(13^), so we can, with the 
current assumption on Q'^^s, consider only a subset where only one of p^i and 
Ps2 is nonzero for each s. 

(ii) The index X„/ with n' = X^s '^s effective charge centers. For X„/, we treat all 
charge centers as distinguishable, so it is computed by the index theorem of the 
previous section with n' < n distinguishable centers. 
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(iii) The combinatoric factor of l/rf^! for each s. This takes into account of the 
reduced volume of the moduh space due to the orbifolding by the permutation 
subgroup S{ds) acting on the reduced n'-center quantum mechanics, but does 
not address the contribution from the submanifolds fixed by S{ds). 

(iv) For each effective particle of charge p/3, with primitive /3 and p > 1, that shows 
up in computation of X„/, one further assigns an effective internal degeneracy 
factor —Q{pl3) = — fi(/3)/p^, in addition to (— i7i)™'(— (^2)''', which reflects the 
fact that m' number of /3i centers and k' number of /32 centers are left as 
individual. 

The last — 0(p/3) = —Q{l3)/p'^ should be compared to the naive {—Q{l3)y de- 
generacy factor that would be the correct factor if we were considering p separable 
particles of charge /3 instead of one particle of charge pP. Finally, the appearance of 
— i7's instead of fi's is natural, since for example a half-hypermultiplet with Q = 1 
acts like fermions, while a vector multiplet with Q = —2 acts like bosons. 

From the quantum mechanics viewpoint, the decomposition (i) clearly has some- 
thing do with the orbifolding action F. Each term in (5.5) arises from a submanifold 
which is fixed by the product of permutation groups of order ps, Yls S{Ps) C. F. For 
each sector, origins of (ii) and (iii) are also evident as coming from a reduced prob- 
lem of n' charge centers and the subsequent volume-reducing action of Yls S{ds) = 
^/ Yls^iPs)- The only part of this formula which is not evident, so far, from the 
moduli quantum mechanics viewpoint is the rational degeneracy factor of (iv). Here, 
we would like to isolate where this comes from, and later derive it directly from the 
moduli dynamics. 

After some careful thinking, it becomes evident that this rational degeneracy 
factor should come from quantum mechanical degrees of freedom normal to the sub- 
manifold fixed by S{p)^s. Let us consider only p > 1 cases and label them ps', since 
otherwise the internal degeneracy factor is fi{f3) as expected. Subgroup S'(ps/)'s per- 
muting these Ps' > 1 charges fixes a submanifold -^n-j](p^/-i) inside A^„. This fixed 
submanifold has a codimension 2 ^(ps' — 1) in A^„, since it is spanned by coincidence 
of Ps' centers, each of which span two directions in A^„. 

Consider the reduced dynamics on the intersection, Ain'=n-J2ip /-i)' ^^^ computa- 
tion of X„/ with all such Ps'Ps' center treated as single particle, respectively. If we start 
with this reduced index problem, impose the statistics, and pretend that the centers 
associated with Ps'Ps' comess with a unit degeneracy we will find a contribution of 
type 



where m'/3i + k' [32 = m[5i + k[32 — Xls'Ps'A' ^^'^ n' = m' + k' + ^^, ds'- Note that 
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we took care to include the volume- reducing effect of Yls S{ds) = ^/ Yl S{Ps') via the 
denominator Yldgl = m'lk'llYlds'l). 

This expression is obtained after ignoring the quantum degrees of freedom that are 
normal to the fixed manifold Ain-J2{p /-i)'^' ^^^ does not agree with MPS formula. 
The latter is 

=i- x„, X i-n,r'i-n,f X n =^ (5.7) 

lls"s- gt Ps' 

so the difference is precisely the rational degeneracy factor of (iv). Clearly it comes 
from quantizing the normal bundles of Ain^-Y^ip /-i)'s inside A^„. On the other hand, 
for all intent and purpose, this part of quantum mechanics is free, since they have 
something to do with many identical particles and has no interaction of type £i, 
except for the statistics issue. *^^ 

This leads us to conclude that the factor, Q{f3)/p'^, should arise from an index of p 
noninteracting identical particles of charge /3, modulo the center of mass part which 
already contributed to Xn'=n-p+i- The relative dynamics of such identical particles 
carry 2(p — 1) bosonic degrees of freedom, 2{p — 1) fermionic degrees of freedom, and 
additionally internal degeneracy of |fi(/3)| for all p particles. In next subsection, we 
will show that precisely such a factor arises from the dynamics of non-interacting and 
identical p particles with the internal degeneracy fl{(3). 

The full MPS formula follows the same set of rules, except that one must in general 
consider an arbitrary set of charges on the -|- side, and all the partitions of the total 
charge 7t in terms of charges of states that exist on + side of the wall. Since the + 
side of spectrum may then include states of charges hf3i + jf32 with h + j > 1, more 
diverse charge centers will appear for the individual index problems on the right hand 
side. As we will discuss later, this can be incorporated by treating all such particles 
on the + side as independent. The only subtlety is when non-primitively charged 
states exist on the + side; this can be remedied by employing the fully general form, 
^(7) — J2p\'y^^il/p)/p'^ 3.S the effective degeneracy factor. We will also see this 
most general emerging from our index computations. 



^^^There is a subtlety, again related to whether threshold bound state of identical charges can 
form. Since we started with the assumption that such nonprimitive state do not exist, it is safe to 
assume this issue does not complicate our problem. Whether or not we can extend this to theories 
with threshold bound states, i.e. supergravity, is an open problem. 
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5.2 Physical Origin of Q,{f3)/p^ from p Non-Interacting Iden- 
tical Particles 

Let us first restrict ourselves to bound states involving several identical dyons of 
charge /3 with —Q~^{P) = ±1.*^^ As in the previous discussion, let us consider the 
bound state of n charge centers, 77- = X^a^-^' "^ '^^ which are /3's. The identical 
nature of the /3 dyons means that the orbifolding group includes the m-th order 
permutation group S{m). We start with the assumption, for simplicity, that kp state 
exists only for |A;| = 1 on the + side, and then come back for the fully general case 
in next subsection. 

Consider the index reduced on the true moduli space Ain as described in the 
previous section, with proper account taken of the bosonic or the fermionic statistics, 

^" (5^^^) - ^^ (5^^^) ^ L^' {{X\i-lY'^e-'^'Vr\X)) dX , (5.8) 

via the orbifolding projection operator Vr- Here tr means the trace over the fermionic 
variables as well as other internal discrete degrees of freedom, and we integrate over 
the bosonic variables X with an appropriate measure. Matching the sign of —Q with 
(—1)''^ value of the component dyon states, as we noted in the case of bound state 
counting in distinguishable centers, this index naturally computes the degeneracy 
0~'s as 

Ind({7A};r)= / tT({X\{-lf'-^e-^'^'Vr\X))dX , (5.9) 

SO we would like to ask whether this reproduces (5.5) and rediscover the rational 
invariant 0. For this, let us concentrate on the permutation group Sim) part of F 
and see how it generates a series of terms, similar to MPS's wall-crossing formula. 

Inside A^„ there are various fixed submanifolds, A^„', of dimension 2(n' — 1). The 
simplest are A^„_p+i, fixed under S'(p) subgroup of S{rn). Note that we use the same 
notation M. for the fixed submanifolds as the full classical moduli manifold A^. This is 
because all of them are of exactly the same type. For example, the manifold Ain-m+i 
would also emerge if we started with a different low energy dynamics involving a single 
center of charge mP in place of m centers of charge /3. Ignoring contributions from 
these fixed manifolds would simply give 

(_l)EA>B(7A,7B>+n-l f^\ X JJ^A (5.10) 

\ I I / A 



"^^^^Because of spin-statistics theorem, there is no irreducible BPS niuhiplet in Z? = 4 A^ = 2 theories 
with n = —1. The half-hypermuhiplet has 1 and vector muhiplet has -2. The assumption here is 
strictly for the illustrative purpose only. 
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due to the volume-reducing action of F when it is acting freely. This is the very first 
term in the MPS formula. Since there are many fixed submanifolds, however, each 
of them will contribute additively on top of this. 

Without loss of generality, let us consider the fixed manifold Ain-p+i associated 
with the partition mf3 = pf3 + f3 + f3 + --- + f3, and label the coordinates along the 
fixed manifold Ain-p+i by X' and those normal to it by Y. Note that among X' are 
the two coincident (or center of mass) coordinates for the p/3 charge center, so we 
can think of F's as the relative position coordinates among these p charge centers; 
therefore there are 2{p — 1) Vs. We then formally write the additive contribution 
from the fixed submanifold Ain-p+i as 

A,xInd({7A'} = {P/9,/9,...};n 
= ApX f tr'(^(r = 0,X'|(-l)2-^3e-^^Vr'|l" = 0,X')) rfX', (5.11) 

where V = r/S{p) is the remaining orbifolding group that acts nontrivially on 
M.n-p+1- Here tr' denotes trace over fermionic and other internal degrees of freedom, 
except those associated with the p identical /3's that are held together at M.n-p+i- 

We factored out the contribution Ap from the normal directions, Y , and the 
superpartners thereof. On the other hand, the second factor is the index of a reduced 
n — p + 1 center problem, modulo the internal degeneracy factor of pfi charge center. 
Other than this, the computation of this latter factor proceeds on equal footing as 

(5.8), 

I tr' {{Y = 0,X'|(-l)2^3e--^Vr'|V' = 0,X')) dX' 

^ (_l)EA'>B'(7^',7B')+n-P f ^n-p+l{{lA'}) \ ^ ""Q ^^, + . . . (5 ^g) 

^ II / A'=2 

SO we may compute the full index recursively. Ap plays the role of the missing 
internal degeneracy factor VLy here, as it computes the effective contribution from 
these p coincident /3's. The ellipsis denotes terms from other fixed submanifold inside 

Mn-p+i etc. 

We will show that Ap = ±l/p^, regardless of precise nature of the (5 particles, 
which also reproduces MPS formula for VL{(5) = ±1 entirely from the dynamics. 
Schematically, this factor can be written as 

~ /"tr f{Y\{-l)^^-^ e-^^^P|r)) dY , (5.13) 
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with Fj_ and i/_|_ defined on F's and the superpartners, again with suitable measure 
for the bosonic integral. The projection operator 

'P = \Y. (t1)''^"^M^ (5.14) 

ensures that we isolate wavefunctions of correct Bose/Fermi statistics. Note that the 
sign in the projection operator is the same as that of — i7(/3). M^ is the (p — 1) 
dimensional representation of tt G S{p)^ common for (p — 1) coordinate doublets Y's 
and for their fermionic partners, -i/^'s. Naturally o'(7r) is odd or even when vr is odd or 
even. 

Since the embedding of A^„_p+i's in A^„ could be very complicated, the exact 
nature of the decomposition is not entirely clear. On the other hand, the initial 
index problem is gapped and allows us to take s — )■ 0. At least in this limit, the 
decomposition makes sense intuitively, and, as we will see shortly it suffices to consider 
an arbitrarily small tubular neighborhood around the fixed manifold M.n-p+i- We 
take the Y directions as a ball i?2(p-i) insides a fiat R^^p^^\ Therefore, we have 



Ap = lim I tr ({Y\ {-lY'^ e^'/^^'^'v\Y)\ dY . 



(5.15) 



Precisely how we cut-off this neighborhood will not matter, as we will see shortly 
that a Gaussian integral of squared width s emerges along Y directions. 

Interestingly, exactly the same kind of object was studied when solving for the 
famous DO bound state problems in the 1990's. The first such computation appeared 
in Ref. [26] on two-body problem and was later expanded to many body case in 
Ref. [27]. Here we will adapt and expand the computations in these works. Since 
there are 2(j9 — 1) real fermions, we will choose a polarization of type {■?/'i, ■i/'l}, so that 
a general wave function | \1') can be expanded as 

I*) = {^{Y) + *^^^(>^)4) + ^^^^'^'H^)4o42) + • • • ) |0) , (5.16) 

where vl>{«i •••*'"} (y) = ^f^X]^^" '^kiX) and {^fc} are complete basis of F-space 
wave functions. Since the Hamiltonian is free and does not mix sectors of different 
fermion numbers, we may evaluate the bosonic and fermionic trace independently. 
The fermionic trace, for each of M^., is given by 

{±lfii^{{-lY-\-lf^^M^) , (5.17) 

where (±1)^ arises from the value of (—1)^"''^ on p individual (5 states. Also the orbital 
part of J^ is always integral in the absence of the minimal coupling contribution. 
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so (— 1)^'^3 acting on the quantum mechanical degrees of freedom becomes purely 
fermionic expression (— 1)*'^^(— 1) '/■''". The sign in front of the latter comes from 
converting the chirality operator to a form involving the fermion number operator 
that counts the creation operators ip'^ . 



Then the contribution from Y direction reads 

neS(p) 



Ap = ^ E (Tin-^r'tr^ ((-1) 



X 

'32(p-l) 



in the limit of s — )■ 0. A crucial observation that allows us to proceed systematically 
is 

tr^((-l)^"M,^ 
= (0|0) - {0\r'M^a'''^l\0) + 1(0|/>'^'M^,,"M.,,VM|0) + ■ ■ ■ 

= det(l-M^) , (5.19) 

and furthermore 

,,/, ,,x fp, Trisa cyclic permutation of order p ,. „^. 

^^^^^-^^) = |o ; otherwise ^^'^O) 

for which it is important to remember that M^ is a p — 1 (rather than p) dimensional 
representation of >S'(p). Since (—1)'^'' = (— 1)^"^ for any cyclic permutation of order 
p, we find 

Ap = ^^ J](Tlf-i(-lf-Met(l-M,0 /(>'|e'^'^'^^'l>')^>' 

7r' 

= ±i y det(l - M^O X — i-- /e-(^-^^'^)'/2^rfr 
p! ^ ^ ' {2ns)P-^ J 

Jl \f det(l - M^O ' ^^'^^^ 

where the sum is now only over the cyclic permutations of order p. There are precisely 
(p— 1)! such permutations and they each contribute 1/p, via the determinant, so the 

result is 

±1 
Ap = — , 5.22 
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as promised. Clearly, we can repeat this when there are several such factors simulta- 
neously, to give, 

^K'}=n^' (5.23) 

reproducing Q{P)/p1 of the MPS formula with Q{(3) = ±1. 

The more general case of fi(/3) = ±(i can be derived similarly. Let us write one 
particle state as 

\^) =-q/\x,ilj)\0;r]) (5.24) 

so that 1] = 1,2, ... ,d labels the internal degeneracy, and p-particles wave function 
(without center of mass degree of freedom) can be written as a sum of terms like 

<ij{^r-k}{vs}^Y)i,l^l^ ■ ■ ■^},)|0;r/i,r/2, ...,Vp), (5.25) 

none of which mixes under the free Hamiltonian. Thanks to this. Just as the fermionic 
part and the bosonic part separately contributed, this internal part also factorizes 
under each vr. Expressing Ap as a sum over the elements of permutation group again, 
we now have an extra factor 

{Vp,---,m\V7r{l),---,V7r{p)) , 

for each permutation it, and the trace over these internal indices. Thus, we arrive at 
a similarly simple form, 

A _ (±1) \r >:^ {vp,---,vi\vnii),---,Vnip)) . . 

^p-^r 2^2^ det(i - M^o • (^-^^^ 

J7l,...,J7p n' ^ 

As before, the sum is over the permutations of cyclic order p. For such tt', the inner 
product vanishes identically unless all ?7j's are equal to one another, and gives unit if all 
are equal. The sum over ?7's thus collapses to a single sum over f] = rji = ri2 = ■ ■ ■ = rjp, 
and 



Ap = (±i)yixi-y^ — - 

p y ) A^ p!^det(l-M^/) 



p2 p2 



(5.27) 



This gives us the only essential ingredient in confirming (5.8), and in fact the fully 
general version thereof. 

For a complete derivation, a recursive argument is needed and we need to con- 
sider possibility of low energy dynamics with charge centers both primitive and non- 
primitive charges simultaneously. This naturally brings us to the most general wall- 
crossing formula, next. 
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5.3 General Wall-Crossing Formula 

Most of what we derived generalizes to cases with arbitrary spectrum on + side, 
without much modification, but here we need to point out one subtlety. Suppose 
the + side of spectrum contains not only a pair of primitively charge states 7 and 7' 
but also states like h'^y + J7' a little more involved, and includes states of composite 
charges such as 7727 or linear combinations with other charges. (One can also have 
states with charges completely unrelated to these but those will not participate in 
the wall-crossing, and therefore irrelevant.) Such a state cannot be considered as a 
bound state of h 7's and j 7"s, since the two are mutually repulsive on the + side. 
Rather, it should be regarded as a completely independent particle of different origin. 
In fact, for SU{2) theory with a single flavor, a monopole 7, a quark 7' and a dyon 
7 + 7' are known to coexist in the central part of the moduli space. 

Let us denote charges of these independent particles as /3„. Since one can form 
bound states of a given total charge 77- on the "— " side from different combinations 
of these + side states. We will label each of these physically distinct combination by 
the upper index inside a parenthesis such that 

1T = Y1 ^v^^v = Y. "^v^^v = J2 "^v^^v = ■■■ ■ (5.28) 

A V V 

In such circumstances, the total degeneracy for 'Jt has to be computed for each of 
such bound state problems and summed over, 

fi-(7T) - 1^+(7t) = 5^^ (K^)}) , (5.29) 

<? 

where each term on the right hand side is computed from the n*^^) = ^ mj center 
quantum mechanics. For each of ^,{{711^ })^s, computation of the previous subsection 
goes through without modiflcation, and will be computed as 

f]({m(^)}) = n\y2^A = Y.<'^M ■ (5-30) 

One important detail to remember is that, even if some of charge /3t,'s might be a 
linear combination of other /3„'s, each of them are physically unrelated independent 
particles: The permutation group is simply T = Ylv^i''^^^ ) ^'^^ each of these index 
problems. 

Combining this with the results of previous subsection, we reproduce MPS formula 
in its most general form. Note that, when we reorganize this formula in terms of the 
index of distinct particles of unit individual degeneracy, 

X„({...,7,...})= / Ch{J^)A{Mn) (5.31) 
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the rational invariants ^(7), multiplying them, will accumulate additive contributions 
of the form, including p = 1 case, 

^^^ (5.32) 

for each 7/p = jS^ that appears in one of the expansion 77- = ^^rnv jS^. Since we 
are summing over all possible such expansions, it implies that 

^W-Z^ (5.33) 

will appear as the effective degeneracy factors that multiply X„'s. p = 1 terms arises 
only when 7 is one of the /3t,'s, while p > 1 terms arise from the orbifold fixed sector 
as in the previous section when 7/p is one of /3„'s. The final expression is 



fi-(7T)-fi+(7T) = 

+ ■■■ (5.34) 



where we wrote a representative form for the partition 77- = X^a'=i 7a' ^^^'^ ^' centers 
and the associated orbifolding group T' permuting among identical elements in {7^/}. 



6 Alternative Derivations for Arbitrary Number 
of Centers 

We have so far seen how Dirac index computes wall-crossing formulae with the statis- 
tics imposed as a projection operator in the index problem. In particular by decom- 
posing the computation into various fixed submanifold, we effectively reduce the 
problem to a collection of index problems for certain sets of distinguishable particles. 
Some of these distinguishable particles originate from real BPS particles, whereas 
some are mathematical construct. In the latter case, say associated with any fixed 
submanifolds under the permutation group S{p), an effective internal degeneracy 
VL[P)/p^ emerges universally. 
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In this section, we repeat this exercise for general cases, starting with a string 
theoretical method of computing orbifold index. While the derivation of previous 
section is very intuitive and compelling, we wish to further support it with this more 
rigorous formulation. To handle this problem, we have to figure out how to define the 
index problem on orbifold singularities. This is a subtle issue in mathematics since 
most index problems are defined on smooth manifolds. On the other hand, in string 
theory the orbifold singularities are common and the string theory on such space is 
well defined. We will use the index for the Dirac-type operator defined in the string 
theory setting and apply it to the problem of our interest. 

In fact in the original paper on the orbifold [45], there appeared the formula of 
the Euler characteristic on orbifolds. The special feature of this formula is that it 
has the non-zero contribution from the twisted sectors and it gives the right Euler- 
characteristic for the blown-up manifold out of the underlying orbifolds when the 
way of the blown-up of the orbifolds is known. This consequently motivates many 
mathematical literatures which try to justify the formulae in a rigorous way [46]. 
For our purpose, the natural object is the Dirac-Ramond index actively studied on 
1980's and we try to read off the index of interest. We are interested in the f/(l) 
equivariant index since we will follow closely the strategy developed by [22]. It turns 
out that the twisted sectors would not contribute and the entire problem effectively 
reduces to the usual Atiyah-Bott-Lefshetz theorem [47] , with string theory serving as 
a computational tool. 

6.1 Dirac-Ramond Operator on Orbifolds 

Now let us turn into the Dirac-Ramond index. Dirac-Ramond index is so called the 
U(l) character- valued index [42] 

/ = lim Tr(-l)^''exp(27rinP)exp(-27rT2H) , (6.1) 

where H is the Hamiltonian for the string worldsheet and P is the momentum of the 
string. Here P plays the role of the f/(l) generator, which is the translation generator 
along the worldsheet. With r = Ti+iT2 and q = e^'^*'^, this index has the contribution 
from supersymmetric states of world-sheet only and*^^ 

I = Tr(-l)^«g^«+^g^o+^" = Tr(-l)^«g--f' , (6.2) 

so that we have the contributions only from left-moving modes. We define H = 
Lq + Lq, P = Lq — Lq and e, e is the zero point contributions from the right and the 



#i8We change the role of q, q in Eq. (6.2) for later notational convenience. 
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left movers respectively. Thus / has the form 

/ = 5]XAg-\ (6.3) 

A 

where X\ is the index on the subspace with momentum A, which is essentially the 
level of the massive string state denoted by m. If we take the lowest A this will give 
the usual index of the Dirac operator. / has the form 

/ = A{R) J2 q'^ChiF, g)Ch{R, t) , (6.4) 

(g,t,m) 

with g representation of the gauge group and t tensor representation of the rotation 
group. The sum is taken over all fields in representation {g^ t) at each mass level m 
of the left mover [44, 48]. 

Now I is closely related to the string partition function. For example, if we 
evaluate the Dirac-Ramond index for the NS fermions with anti-periodic boundary 
conditions along worldsheet time direction in the heterotic string theory, we have [43] , 






where Xq, ip'^Q denote the bosonic and fermionic zero modes. When we set F = i? = 
we have the typical form of the string partition function. In the above we defined 

R,^u = -Rf^,yXaixo)^p\tp% , 

Fab = —F^^T^^r-or-o , (6-6) 

and regard the fermion zero modes as the differential forms in the usual way. For our 
purpose, the gauge bundle is Abelian, so we simply replace F here by J-' at the end 
of the computation. 

If we take the lowest q this picks up the basic representation and / gives 



dx!^dtPtoA{M)Ch{F) , (6.7) 



M 



where we use [44] 



y^x''Ch{F,[k]) = det{l + xexp^, 

y"x''Ch(F,{k)) = det{l - xexpl^)-\ {Q.i 



41 



with [k], {k) being the antisymmetric and the symmetric products of the basic rep- 
resentation. Thus in this way we obtain the formula for the Dirac index from the 
Dirac-Ramond index. 

Now this formahsm can be easily generalized to the orbifold cases. The twisted 
index of the Dirac-Ramond operator is given by [49] 

Iig^,g2■,r) = TTg,[f'+^q^'^'g^i-lf-] (6.9) 

where the gi,g2 are the twists along a2, Ci directions, respectively. Note that we define 
gi twists along temporal direction. Here we label the points of the world sheet torus 
by the complex quantity ci + ro'2 with cti, o"2 having the periodicity 27r. Consistency 
of the boundary condition requires that gi, g2 commute. In this case the bosonic and 
the fermionic zero modes exist only along Mg^g^, which is the space inert under gi 
and g2- If we decompose the curvature along and the normal directions of dimension 
dn to Mg^g^ whose dimension is d, 

^Rm = ^i?.H,^_o^io = ^ _,,^ J ' ^ = ^ ■ ■ ■ 2 ' 

while gauge bundle of rank d^ has the field strength components 

(6.10) 

where i,j, k, I denotes the indices along M^^^j directions. 

Let's also define the action of gi, (72 on a real bundle V. V could be either tangent 
bundle TM or the gauge bundle E. Under the action of gi,g2, V restricted to Mg-^g^ 
formally decomposes as 

y\M,,,, = y^'^ ® i®rV^''^)- (6-11) 

The second term represents a formal sum of 2-dimensional real bundles on which 
91,92 act as rotation by the angles 2nar,27if3r respectively. Then, 

1(91, 92; r) = / dx^d^pl 



det'{d + i^ 



. H — —^ . (6-12) 
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where det' denotes the determinants without the zero modes. Here we denote deta,i3{d+ 
v) denotes the determinants of (9 + z/ of complex modes ^(ci, 0-2) with the boundary 
conditions, 






(6.13) 



And 



deta,i3{d) = d 



a 



r 



2_fl , 1 






JJ(1 - g"-^e^"^")(l - g"+/3-ie-2"^") , (6.14) 



n=l 



with ^ a,/? < 1. Finally 

deta,p{d + v) = e~'^^''d 
If we evaluate one such component 

deta,ii{d + ipr)det_a-i3{d — ipr) 



a + u 

(3 



(0|r) . 



(6.15) 



e 2 



(2a/3-a-/3) 



(^-hy 



q 



n(i 



n=l 



^n-/3g2^iag-2^p,^)(^^_^^ 



.n+l3—l—2nia2npr 



for /3 ^ while for /3 = 



Jdeta,o{d + ipr)det-afi{d - ipr) 

00 
e-x"g-^ JJ(1 _ gV™e-2^^'-)(l - g"e-2^*"e2^'''-) 



n=l 



X . /Q p2iTiap~2ivpr\(Y p~2-nap2ivpr\ 

Note that the zero mode contribution is given by 



^J{1 - e2'ri"e-2'^/''-)(l - e-'^'^^^'e^'^P-) = ±2isinh7r{pr - ia) 



(6.16) 



(6.17) 



where ± sign can be chosen so that it reduces to \2simia\ with pr = 0. Given this 
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machinery, one can evaluate the twisted index of the orbifold 

1 



/(fifi,0:r)= / dxndth\AM„ — i 

nil nr=i(l + g"-le2™-e-2^^'-)(l + g"-le-2^*"'-e2^'''-) 1 






-27riapp27Tcjp) ±2isinh7i{ujp - iap 



d„ 



rfx^dV' ^ — 4-. 7^-^^ TJTT TT r^ r 6.18 

where |^, ^, ^ have the matrices with skew eigenvalues are co^, cOp — iap, pr — iotr , 
respectively. For V'^^'' directions, ar = is understood. In particular we are interested 
in the basic representation and we will read off the term of the lowest order in g, which 
is g2 in this case, 

dn 

f dxld^pi,AiM,Jexp'-f- n . ] — = Tr(-l)^^i , (6.19) 

JMg^ Svr f^ ±2tstnh-K{ujp - tap) 

which is nothing but the Atiyah-Bott fixed point formula. In writing down the action 
we discarded the phase factor depending on a and the zero point energy contribution 
to the fractional factor of q to match the usual index result. Note that sink factor 
gives the gi action on the normal component of the curvature. 

Now for general twisted sectors 



-'9192 n= 



\ det{l - q^exp'£) 

ii. 1 1 



JMa^a. f-1 detil - q^'exp^) 

_ (6.20) 

TT i. TT°° (1 — (7"'-/5pe^'^*"Pe~^'^^p)(l — (7"'+/3p"le-27riapg27rajp'\ 

Here F is the field strength along V'^'^^ while d^ denotes the dimension of the vector 
bundle excluding V'^'^\ The peculiarity of this expression is that there will be contri- 
butions only from the states invariant under the orbifold action. One can easily check 
that the projection operators — ^^ ^^ ^^ — for an element gi of order m acting on 
the twisted sector by g2 gives the vanishing contribution to the twisted index for any 
state not invariant under the orbifold action. Thus one of the nontrivial contribution 
comes from 

/ MMg,g,)qhxp'f , (6.21) 
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while more complicated combinations of the tensor representations of the normal bun- 
dle of the tangent bundle and V'^^^ are possible. The above term Eq. (6.21) comes 
from the twisted states localized along Mg-^g^ and singlet under the rotation group of 
the normal directions to Mg^g^. If we consider more general states transforming non- 
trivially under the transverse rotation group, they have more complicated expression, 
which can be read off from Eq. (6.20). 



6.2 Rules for Twisted Sectors 

Since we could in principle have the nontrivial contributions from the twisted states, 
we have to decide if we have to worry about contributions from the lowest twisted 
states when we are trying to read off the index for the field theory limit. We already 
know that for Euler characteristic we have to include the contribution from the twisted 
sectors. We also expect this is the case for the signature of a four manifold which 
counts the difference in the number of self-dual and the anti-self-dual 2-forms since 
two forms could arise in the twisted sectors. 

To facilitate the discussion, let us consider a specific compactification of string 
theory and figure out the plausible rules for each index problem. Consider Type 
IIB compactified on K3. This is a chiral N=(2,0) theory and the matter contents 
are completely fixed by anomaly constraints. Thus for whole moduli space of K3, 
the spectrum is the same and this applies for K3 orbifold such as T'^ / Z2-,T'^ / Z^. 
The matter content consists of 1 supergravity multiplets and 21 tensor multiplets. 
Supergravity multiplet consists of graviton and 5 self-dual two-index tensors and 
gravitinos. One tensor multiplet consists of 5 scalars and one anti-self-dual tensor 
and chiral fermions. The tensors arise from two massless states -B* with i = 1,2 
of Type IIB string theory, and compactify self-dual 4-form -B^^;^„ on two forms on 
K3. On K3 we have 3 self-dual 2-forms and 19 anti-self-dual 2-forms so that we 
get 3 self-dual tensors and 19 anti-self-dual tensors out of B'^^-^ . If we realize K3 
as T^/Z2 orbifold where Z2 acts as Zi — )■ —Zi, Z2 — )■ —-22, for a complex coordinate 
Zi of the torus, we have 16 fixed points and, upon the blowup, each of the fixed 
point is replaced by Eguchi-Hanson space, which supplies one anti-self-dual 2-form. 
Hence on each of the twisted sector we obtain one tensor multiplet. Now from this 
consideration, it's obvious that we have to include the twisted sector contribution for 
the computation of the signature of K3. This is given by difference of the number of 
self-dual tensors and anti-self-dual tensors. 

Now consider a Dirac spinor of positive chirality spinor ri G 6*^3 on K3. If we 
have such spinors, one obtains the 6-d gravitino from 10-d gravitino 

,M 



^'^' ^ ^M ® r/ (6.22) 
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Thus if we have such spinors on K3, one can obtain gravitinos on 6- dimensions as 
many as those. One can ask if such spinors could arise from the twisted sectors. But 
we already saw that on the twisted sectors, only tensor multiplets arise so that we 
do not have gravitinos from the twisted sectors. Had we obtain the gravitinos from 
the twisted sectors, we will also have multiple gravitons by supersymmetry, which is 
quite weird! Thus spinors of K3 can come only from untwisted sectors from the Z2 
invariant projection of spinors on T^. These spinors are coming from Ramond sectors 
of the right moving sectors*^^ and we can see why the spinors of K3 could not arise 
from the twisted sectors. In the untwisted sector 8^ of 10-d postive chiral spinor can 
be represented as Ramond vacuum 

|±,±,±,±>, (6.23) 

in the light cone gauge with total number of + sign being even. Here four ±'s 
denote the eigenvalue of exp{i7TJi) where Jj are 4 Cartans of 5*0(8). Under 5*0(4)1 x 
50(4)2 = SU{2)li X 5f/(2)ra x SU{2)l2 x SU{2)r2 with 50(4)2 action on the 
transverse K3 directions 

8^ ^(2, 1,2, 1)0(1, 2, 1,2), (6.24) 

and under the Z2 orbifold action only (2, 1,2, 1) is invariant. *^° Combined with 8^ 
of the right moving modes transforming as a vector under 5*0(8), 

8,^(2,2,1,1)0(1,1,2,2), (6.25) 

where (2, 2, 1, 1) is invariant under the Z2 action. The combination 

(2, 1, 2, 1) ® (2, 2, 1, 1) © (1, 2, 1, 2) ® (1, 1, 2, 2) (6.26) 

give rise to the gravitinos of 6-d after Z2 projection (plus fermions of spin |). Now 
consider the twisted sectors. The fermions arise from Ramond vacuum of right- 
moving sectors possibly combining with the NS sector states. Important thing is 
that Ramond vacuum of the twisted sectors is given by 

|±,±> (6.27) 

with the total number of + signs even. Note that this state transforms nontrivially 
only under 5*0(4)1 and is a singlet under 5*0(4)2. This could not be a spinor on K3, 
which transforms nontrivially under 5*0(4)2. This argument goes through exactly for 
the other K3 orbifolds. 



"^^^Spinors could arise from the left moving sectors for Type II string theory but the arguments 
goes parallel to the right-moving modes. 

'^'^^Z2 action on Ramond vacuum can be written as expiir^J^ — J4)) where J3, J4 are the Cartan 
elements on SO{4)2- 
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For the situation we are interested in, we have to deal with a spinor in S'^ ® E 
with a suitable vector bundle E. The associated Dirac index is naturally defined 
in the heterotic theory. Since the orbifold action acting on the right-movers cannot 
generate S~^ one cannot obtain the desired spinors from the twisted sectors. 

This argument goes through Calabi-Yau 3-fold case as well since again the Ramond 
vacuum of the twisted sector is of the form 

|± > (6.28) 

and this is a singlet under 5*0(6) of the transverse direction of Calabi-Yau 3-fold 
while spinors on Calabi-Yau 3- fold transform nontrivially under 5*0(6). Hence in 
the spacetime supersymmetric theory, we do not have to consider the contribution 
from the twisted sectors for the Dirac index. This is also true for nonsupersymmetric 
orbifold, since these states in the right-moving modes are the excitation with the 
tensor representations of the transverse rotation group around the Ramond vacuum 
of the type Eq. (6.28). This could not match the spinor representation of the spinors. 

Now let's apply this formalism and carry out the Dirac index computation for the 
simplest orbifold T^/Z2. One can use the formula Eq. (6.19). From the projection 

^Tr ((-1)^(1 + «)) (6.29) 

with a being Z2 action on T^ we find that the first term |Tr(— 1)^ is vanishing since 
this counts the Dirac index on T*. From the |Tr(— 1)^(7 one obtains 

-■16-- — ^—- = 2, (6.30) 

2 (2smf)2 ' ^ ^ 

where 16 comes from the total number of the fixed points and 2sm| factor comes 
from the Z2 action on the normal bundle over the fixed point. Indeed we obtain the 
right Dirac index on K3. Similar computation can be done for the other toroidal 
orbifold of K3 and gives rise to the same result. 



6.3 Orbifold Index and the Equivariant Generalization 

Let us come back to our problems of counting BPS bound states of charge centers 
and see how the Bose/Fermi orbifolding is refiected on the equivariant index. Con- 
sider the simplest case where we started with a pair of identical particles and other 
distinguishable ones. The orbifolding group would be then 52 = ^2- Here let us 
assume that —Q~^ = ilj and then the expected answer, from either MPS formula or 
the derivation in the previous section, is 

^(_l)(27i,72)^(27i + 72) = ^/(7i,7i,72)±/(27i,72)^ • (6.31) 
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On the other hand, we may also read off the index from the ^(1) character- valued 
index, with the projection operator 

1 ± a , ^ 

^- , (6.32) 

with a^ = 1. A new ingredient here is the sign in front of a, which chooses whether 
one retains bosonic or fermionic wavef unctions. 

When we apply the formula derived earlier, we have to know how the bundle of 
interest would transform under Sn actions in general. Note that we are dealing with 
rank one bundle given by 

-^ = -^ E '<^bc ""T'^iZ QAB dx^^ A dx^K (6.33) 

2^;^ \xa-xb\' 

One can check that the vector bundle is invariant under the exchange of any two 
xa, xb and hence is invariant under the action of Sn- Reading off individual terms 
from the formula in the last subsection, the result for 5*2 = Z2 is 

I [ dx^,dij^_,A{M)Ch{F)±l [ dxldij\A{M2)Ch{F) ^. .\ , (6.34) 
^ Jm ^ JM2 Zismrnxuji 

where M2 is the manifold fixed by Z2 action. The first term produces the first 
term in Eq. (6.31), so let us take a look at the second term, oo'^ = Ui — iai with 
Ui is the eigenvalue of the curvature tensor along the normal bundle and 27rQ;i is 
the rotation angle of the Z2 action, which is vr. Here we assume that the internal 
degeneracies are taken into account as explained in the main text. For the 3 body 
case with 2 identical particles ui does not contribute for the dimensional reason 
since 2isinh'Ku[ = 2cosh'Kooi gives only even powers of ui and we have the 2-d fixed 
manifold with 2-d normal bundle, sinh factor is due to the Z2 action along the normal 
bundle, which gives rise to 2sln^ ~ V^- Combined with 1/2 factor in the projection 
operator this gives the wanted ±1/2^ factor at the second term of Eq. (6.31). 

For 5*3 orbifold and beyond, it is more convenient to consider the U{1) equivariant 
index and use the fixed point theorem where U{1) is generated by the third component 
of angular momentum operator J3. For general action of ^f,*^^ 

Tri-l)V^' = I, <(.irfV'' „-4(M,. ^)CMF. .) n ^2.s,nhU^'^;:) ' ''■'" 



'^^^"'^Here we have j/^tt,/;? -y^jijig jj^ i\yq previous section we use y"^"^^ in the index definition. This 
facihtates the comparison with the corresponding formula appearing in [23] . 
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where u'^ = Up — ia^ with ojp is the eigenvalue of the curvature tensor along the 
normal bundle and 2'Kap is the rotation angle and e^ = y. Following [23], we change 
the expressions into equivariant characteristic classes accordingly and L denotes the 
action of J3 on T*^^'°^(M„). Now consider S2 = Z2 case. This equivariant index can 
be evaluated by evaluating the fixed points under J3. The fixed points are simply 
north and south poles of the sphere, L has the eigenvalue ±1 and a;^ = assuming 
we are dealing with isolated points. Thus, 

-^ 7T--Z77. ^ = TT— 17 = ..,. ..-9 ■ (6-36) 



i2sinh{vL + Tiu'i) 2ism/i(z/ — y) 2coshv y'^ — y "^ 

Note that the expression is invariant under u — ?> —u so that we don't have to worry 
about the sign of L = ±1. This persists for all Sn case as well. 

Now consider S3 case. S^ projection has the form 

^{l+g + g^){l±a) , (6.37) 

where g = (123) meaning xi — )■ 0:2, a;2 — )■ 3:3, X3 — )■ Xi. Again, we inserted a sign 
± in front of the order 2 permutation a to pick out either bosonic or fermionic 
wavefunctions. Then g'^ = (132) and a = (12), meaning Xi — > a;2,a;2 — )■ Xi. One 
can check that ga = (13), g'^a = (23). For the bulk contribution without g or a, the 
projection operator simply gives 1/3!. Now consider 

±^{a + ga + g^a) . (6.38) 

Each of the three elements give rise to Z2 fixed manifold which we already worked 
out. We have the contribution 

3! y^ — y ^ ly^ — y '^ 

Now turn to hig + g"^)- These two elements are order 3. For each element we evaluate 
the equivariant index again by going to a fixed submanifold. The A genus combined 
with Chern characters give rise to the usual angular momentum factor. The relevant 
factor for the orbifolding effect is 

-1 



y-y 



2sinh{u + f) 2sinh{u - f) y^ - y-^ 
Thus, for the refined case 

-(-l)'^^^'^^^fire/(37i + 72) = ^X,e/(7l,7l,7l,72) 



(6.40) 



±Iref{2ji,-fi,j2)-— 3^+Xre/(37i,72)-^ ^ , (6.41) 
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including the factors of the internal degeneracies. Here we are considering the moduli 
space of three particle configurations invariant under 6*3. Normal directions mean the 
deformations away from the fixed configurations in Z3 invariant way. The action is 
rotation by ±27r/3. Taking y = 1 limit, we find 

-(-l)'^^^'^^^fi(37i + 72) = ^/(7i,7i,7i,72)±/(27i,7i,72)^ + /(37i,72)^(6.42) 

reproducing the ordinary index formula of the previous section. 

It is now obvious how to generalize for general permutation group. Again we are 
assuming —Q'^ = ±1. If we consider the various manifolds fixed under the elements 
of order m the contribution from the normal directions in the index formula gives 

1 
2sinh(u + ^)2sinh(u - ^)2sinh(u + ^)2sinh(u - 2^) ■ ■ ■ 

= -^^^^ , (6.43) 

yin y^™- 

which becomes 1/m in the y — )■ 1 limit. Consider m(5i + P2 for general m. One has to 
figure out the combinatoric factors appearing in the projection of Sm group elements. 
Suppose that m can be written as 

m = niirii + ■ ■ ■ + niini. (6.44) 

The number of ways of the partition of m elements into rii of nii elements is given by 



I 
m 



J. J. ri fm I 



nAm. 

i=l 



For example if we partition 4 into 2 of 2 elements, we have 3 possibilities, i.e., 

12|34, 13|24, 14|23. (6.46) 

For each particular partition of m elements, one has nj=i(^^« ~ 1)' permutation ele- 
ments. For example if we consider the particular partition of 6 into 123|456, from 123 
we have two permutation elements (123), (132) and from 456 we have (456), (465) so 
that one can generate 4 permutation elements out of this particular partition. Hence 
the total number of the permutation elements arising from the partition eq.(6.44) is 






n ^- (6-47) 



nAm. 

1=1 
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Combined with ^ in the projection operator of Sm group, and the contribution from 
the normal directions in the index theorem, we obtain the factor 

ri^ri(- i"C )" (6.48) 

Thus for simple case of — f2+ = ±1 we have 

= E n ^^-/(mi/3i, mi/3i, ■ ■ ■ m,A, m^) {[i— f ~ ^ \_ ■ -n,r{-n, 

^-^ -'■-'■ n,- ! -'--'■ m,- y™-' — y m, 

where Ire/ has rij of mj/3j factors. The sign of —Qi can be taken account by the 
judicious choice of the signs in the projection operator of 5*^ as was done for 5*2 and 
5*3 cases. Thus we obtain the expected answer. 



7 Summary and Comments 

In this note we showed how n generic BPS dyons of Seiberg-Witten theory interact 
with one another, and how the relevant low energy dynamics with A/" = 4 supersym- 
metry can be derived in the vicinity of a wall of marginal stability. The resulting 
quantum mechanics is specified by three classes of quantities: kinetic term, poten- 
tials, and minimal couplings. The latter two turn out to be constrained to each other 
by supersymmetry and can be derived exactly, and are universal, in that the general 
structure is applicable to BPS black holes as well. The kinetic term may differ, but 
for counting non-threshold bound states via index theorem, we only need the asymp- 
totic form of the kinetic terms, which fixes effectively the entire Lagrangian. Thanks 
to the universal form, this Lagrangian can also be used to compute non-threshold 
bound states of BPS black holes as well as those of Seiberg-Witten dyons. 

We showed how the usual truncation (in the previous BPS black hole studies) 
down to zero locus, A^„, of potentials is misleading because the massgaps along the 
classically massive direction are always the same as the quantum massgaps along 
Ain, due to the latter's finite size. Instead, one must sacrifice A/" = 4 supersymmetry, 
in favor of an index-preserving J\f = 1 deformation, in order to reduce the problem 
to a nonlinear sigma model on A^„- 

This gives a definite prescription, hitherto unknown, on how to handle the fermionic 
superpartners, and the final form of the index is that of a Dirac operator on A^„ with 
an Abelian gauge field J-" determined unambiguously by the minimal couplings among 
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dyons/black holes. Along with n — 1 radial, classical massive directions, 2(n — 1) 
fermionic partners become decoupled from the problem, leaving behind a supersym- 
metric quantum mechanics on A^„ with real supersymmetry. (Three bosonic and 
four fermionic variables decouple also, playing the role of the center of mass degrees 
of freedom.) This shows rigorously why the Dirac index is the relevant one, as was 
anticipated by de Boer et.al. [25]. 

Since typical wall-crossing problem involves only two linearly independent charge 
vectors, and thus bound states of many identical BPS states, statistics is of major 
importance. We address this directly for the index problem by inserting the relevant 
projection operator Vr, and expanding the index to a series involving various fixed 
submanifolds. Each such contribution consists of two multiplicative factors: one is 
usual Dirac index on the fixed submanifold and the other is contribution from the 
normal direction. The latter turns out to be universal and generates a numerical 
factor ~ 1/p^ for each p coincident and identical particles, times the intrinsic de- 
generacy of the particle in question. This eventually lead to the rational invariants, 
^(7) = Sp|7 ^{1 1 V) I V^ ^ cis the effective degeneracy factor, as was also noted by Man- 
schot et.al. [22]. In the end, we have derived the general wall-crossing formula, from 
the viewpoint of spatially loose BPS bound states by starting from Seiberg-Witten 
theory, ab initio. "^^^ 

Another unexpected bonus, which also shows why the low energy quantum me- 
chanics is never really 2(n— 1) dimensional but must be addressed with all 3n position 
coordinates, was a clarification of various index quantities in literature. The field the- 
ory index is most generally computed by the protected spin character, 

Tr ((-l)2-^^^y2J3+2/3) , 
We showed how this quantity maps to the quantum mechanical index 

Tr ((-l)2^3^2J3+2/3^) ^ 

with the latter's 5*0(4) = SU{2)l x SU(2)r invariance. This is then further reduced, 
mathematically, to the usual equivariant index defined over Ain as 

Tr((0| (-1)2^3 |0),ea.,l/''^^) , 

where jTs = J3 + /s is naturally picked out because in the reduction process, which 
is a purely mathematical operation, we lose part of the 5*0(4) = SU{2)l x SU{2)fi 



"^^^The only unresolved issue here is the angular momentum content of the bound states, which is 
in particular needed when one translates the quantum mechanics index to the second helicity trace 
of D = 4 and N = 2 field theory. For this, we followed the well-known assignment, which has been 
tested in many explicit examples and widely believed to be correct. 
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symmetry. At the end of the reduction process, only SU{2)j survives and rotate 
bosons and fermions equally, unlike the physical rotation symmetry SU{2)l which 
sees that fermions arise from D = 4 spinors. Also we saw how the chirality operator 
(^_X)2J3 ig related to the canonical one (— l)-^-^^ as 

Tr ((0|(-l)2-^«|0),ea.,y''^0 = (-l)^-<-<^-^->+"-^ X Tr {{-if-^-y'^^) , 

giving us by-now familiar sign in wall-crossing formulae. This formula is valid before 
statistics imposed. Imposing Bose/Fermi statistics introduces a further complication, 
which are addressed in much detail in section 5. 

An important question that remains unaddressed satisfactorily here, or in any 
other existing literature, is whether and under what circumstances threshold bound 
state appear. Primary examples would be states that are bound state of more than 
two of the same BPS states. In all of D = 4 A^ = 2 field theory examples where explicit 
states were constructed, no such states ever appeared as far as we know, yet BPS 
black holes can be easily found to have an integer multiple of a given charge. This tells 
us that, for addressing such problems, details of the kinetic term in the low energy 
dynamics are important and perhaps one should deal with 3n dimensional low energy 
dynamics directly, instead of reducing the problem to A^„; A^„ would be no longer 
compact and thus less useful, since classical ground states would include configuration 
with arbitrary small separations and those with arbitrary large separations among 
some charge centers. We further expect that the kinetic terms for dyons and for 
black holes, respectively, will be substantially different at small mutual separations, 
accounting the difference. For dyons, this poses additional difficulties since we must 
face the fact that dyons in field theory are actually non-Abelian objects. For black 
holes, horizon effectively acts as short-distance cut-off and the configuration remain 
Abelian at any separation, and the question of kinetic terms at small separation is 
better posed. We hope to return to this problem in near future. 
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Appendix 



A A/^ = 4 Superfield Formalism 

According to [38] and [37], one can introduce A/" = 4 superfield by 

^^p = iD^Df, + D^D^)V (A.l) 

with V being the D = 4 N = 1 vector supermultiplet. Following [38], we define*^^ 

n - A_lfl ^ 

5" = -^--ea-^ (A.2) 

89^ 2 dt ^ ' 

so that 

{D^,D^} = -t6f^^ (A.3) 

and the indices are raised and lowered by e^/j = — e"'^. Alternatively, ^a/3 is uniquely 
determined by the conditions 

$«/3 = $"^ , 

Da^,3j + £'/3$7" + Dy^l^a = . (A.4) 

It is convenient to define real superfields $ by 

i_i-l- 1 
= a:a + -Oo-aX - 7pC(^aO + -OaaOF + -eabcXb9ac6 

8 8 lo 

This defines A/" = 4 superfield for 3 bosonic and 4 fermionic degrees of freedom. The 
supersymmetric transformation is given by 

X i _ i _ 

oxa = 2^(^aX - 2^^^*^ 

6X = -eCTaXa + -F 

5F = ejc-ex. (A.6) 



#23rp]^jg jg ^ different convention from N=l so we may have to cliange tlie normalization of the 
fields accordingly to conform toN=l convention used in the text. 
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The kinetic term is given by 

1 



Co = -J d'OLi^a). (A.7) 

With / = dadaL{xh) this can be written componentwise*^^ 

^0 = f{X){\xJ-'-xx + lxX+lF') 

1 1 • 1 

-^FxdaXdaf + -^eahcXb{x(^cX)daf - -^XXXXdlf ■ (A.8) 

Now the generahzation to n superfield case is straightforward. One introduces n 
superfield <l>^, A = 1, ■ ■ -n. The kinetic part is given by 

Co = ljd'0mt) (A.9) 

This action is rather comphcated in component form. 

Note that in this simple A/" = 4 superspace formulation, potential terms are not 
obvious, *^^ which is reasonable since the superfields $'s came from the gauge vector 
multiplet of D = 4 by dimensional reduction. Furthermore, as we emphasized in 
the main text, J\f = 1 superspace description is more convenient since the evaluation 
of the index actually relies only on A/" = 1 supersymmetry. For this, we split the 
superfields $ to 3n bosonic superfields and n fermionic superfields as 

<l>^ = x'X-ie^'X,A'^ = iX^ + ieb'^, A = l,---n, a = l,2,3. (A.IO) 

where the complex fermions x and x are decomposed into 3+1 real fermions ip'^^s and 
A's, and the auxiliary field F is now redefined as 6 = F/2. With these superfields, 
the general form of the kinetic Lagrangian Cq was worked out in Ref. [28] 

Co = jde (l^f^D$f $f - UabA^D^^ - ^/^^$f A^ + ytBcD<^aD<^!D^c 

+^<Bc^<^f ^<^f A^ + ^m^Bc<^a A^A^ + ^/^BcA^A^A^] (A.ll) 



"^^^^A simple version of this with constant / would be the dimensional reduction of supersymmetric 
QED [50]. 

'^^'''Nevertheless, Ref. [51] offers a superconformal example with potential terms, written in an 
A/" = 4 off-shell form. 
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with D = -^ — i9^, where 



9ab = {S^A + e^'''^fe)dld^BL 

hAB = Sabd'^d^L 

fa _ anh oc r 

JAB — ^bc '^A'-'b^ 

1 

r^abc _ ^pqh a b cp^a oh Qc r 

^ABC ~ 2 '^pl Qm'^hn^A^B^C'^ 



1 
„ab _ (f:Pqn a b 

"-ABC ~ 2 ^ P'' 9"^ 



e an'^nife - e™ „(5ia)(9^(9^(9^L 



™a _ j aoaobocT 

'"'ABC ~ 2 ™" i' '^a'^b'^c-'^ 
Use = ie,,,919|9^L. (A. 12) 

Note that, for the asymptotic form L{x) = 'Y1,a''^a{,xaY /'^^ all but the first two, g 
and /i, vanish identically. 



B Reduction to Nonlinear Sigma Model on M. 



n 



With n centers, one starts with 3(n — 1) bosonic coordinates and 4(n — 1) fermionic 
ones, after the free center of mass part is removed from the dynamics. It is convenient 
to work with a coordinate system where n — 1 of them equal to independent linear 
combinations of /Ca's. In a slight abuse of notation we will denote these again by /C"^, 
now with A = l,...,n — 1; although there are n /C's, only n — 1 of them are linearly 
independent. Thus, we split the relative part of r^" and -^^^ as Z^ = (/C^,|/^) and 
ip^ = {ip^, "ip^), with M = 1, . . . , 3(n — 1) and // = n, . . . , 3{n — 1). Along the same 
spirit, we also denote by A"^, n — 1 linearly independent combinations of A's that 
belong to the relative part of the low energy dynamics. What do we mean by ■ip'^7 
We wish to preserve at least one supersymmetry, say (^4, and naturally ip'^ is the 
superpartner of Z^ , 



^'' = ^Il;^^^^ (B-i) 



dz^ 

and the kinetic term of ip^ includes two factors of dr^"" /dZ^ . 

As argued in section 4, it suffices to consider the dynamics with flat metric, which 
after taking out the center of mass part becomes 

QAaBb = ITT-AB ^ab , 

where itiab is the (n — 1) x (n — 1) reduced mass matrix. Expressing this in the 
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curved coordinate system, Z^, 



E 



QMN — 7 ^f^AB 



' dZ^ dZ 



N 



we find that partial derivatives of metric coefficients Qmn are nontrivial. In contrast, 
notliing much happens to A's, other than one of them being taken out as the center 
of mass part, so their metric is the same reduced mass matrix, 

hAB = mAB , 

and is constant. Thus, no coordinate-dependent transformations are needed for A's. 
The deformed Lagrangian with flat kinetic term reads in this coordinate, 

C = ^gMN{Z)Z''z''-^e{m-yiCAiZ)ICB{Z)-W{Z)MZ'' 

+ '- (7Miv(^)^^^^ - '- dL9MNiZ)Z''^p^^p^ + '- uiabX^X'' 

+ lidBlCAil^^'X^ + idMy\^N{Z)^'^^'' (B.2) 

where the crucial middle term in the second line follows from Eq. (B.l) and the 
anticommuting nature of fermions. We also used O^ICa = 0. 

Since we anticipate that /C directions will decouple as .^ — ?> oo, we split the metric 
as 

[9a 



_ Hab Ca^i 

MN\ — \ ^T r^ 

and the likewise for its inverse 



.MNi_i {H-CG~^C^)-^ -{H - CG^^C^y^CG-^ 

i9 \-\ -G-^C^{H-CG-'C^)-' G-^ + G~^C^{H-CG'^C^)~^C^G-^ 

Ignoring W and fermion contributions to the conjugate momentum for now for sim- 
plicity, the bosonic part of Hamiltonian will then looks something like 

1 1 

= \{G-Tp,p. + \g^''PAPB + --- (B.3) 

where 

Pa=Pa+{H- CG-'C^)ac9''^P, =Pa- iCG-')/p^ 

Pa's have the standard canonical commutator with /C's, so it is clear that, together 
with ~ ^^JC^ terms, they form very heavy harmonic oscillators of frequency ~ ^, settle 
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down to its ground state sector, and decouple from ground state counting. This leaves 
behind 

nc^^iG~Tp,p. + --- (B.4) 

Denoting the canonical conjugate of j9^ in this reduced dynamics again by y^,*^^ the 
corresponding Lagrangian would be 



^ — c) ^A"^ 



rr + ■■■ (B.5) 

K=0 



This makes clear that we could have done the same more simply by imposing /C = 
at the level of Lagrangian. 

Procedure leading up to (B.4) can be repeated in the presence of Ws, which 
simply shift the conjugate momenta in the Hamiltonian, and it is clear that only 
W^'s will survive. We should ask whether this is consistent, since after all (iW's are 
Dirac quantized magnetic fields, and removing some part of the gauge connection 
could make the remainder ill-defined. However, we have 

dW = deWAdlC^dlC^ + {O^Wa - <9aW^) dy^'dlC'^ + d^W^dy'^dy^' 

and the pull-back onto A^„ is simply 

M'.idW) = d.W^dy'dy^ (B.6) 

The pull-back of a well-defined bundle to a smoothly embedded submanifold is still a 
well-defined bundle, so the reduced gauge connection W^iJC = 0) is consistent. Thus, 
the bosonic part of the action reduces to 



2 ^ 



rr - >v^ 

/c=o 



r + --- (B.7) 



leaving us with the question of how to reduce fermion sector. 

The fermions enter the Hamiltonian in two places. One is as bilinear connection 
term added to the conjugate momenta, and the other is an additive contribution of 
the form 



^^^Generally K. will mix in the definition of this new y^ coordinates, to reflect the shift of the 
conjugate momenta, but this becomes irrelevant because dynamics forces /C = 0. Therefore, the 
same old y coordinates can be used here. 
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with the canonical anticommutator among ■j/'^'s equal to g^^^ . To disentangle heavy 
ip^ from light ip^, we shift the light fermions as 

such that 

Let us categorize these fermion bilinears into three difference pieces, 
-td^W.^'^r - iSa^^^^V + [-t^dslCA^P^'X^ + ■■■] ■ 

Terms in the last bracket involve only ip"^ and A'^'s with eigenvalues ~ ^, so these 
will decouple from the low energy spectrum. The potential mixing between heavy 
and light modes are in 

Sa, = OaW, - d,WA + {CG-')/{d,W^ - d^W,) . 

For heavy sector, this is of course a minor perturbation and ignorable as ^ — )■ 0. 
For light sector, things look less innocent since the size of this operator is itself not 
negligible. However, the heavy fermion enters this operator linearly, and always will 
connect excited states and ground states of heavy fermion sector. This forces the 
energy eigenvalue differences {En — Ek) in the denominator of the perturbation series 
to be of order ~ ^, such that the perturbation is suppressed by powers of ~ S/C,- In 
the end, again, the net effect is to turn off the heavy modes ^'^ and A"^ completely, 
leaving behind 

only, where we will call this light fermion again as ip^^s. The simplest way to under- 
stand this is to recall that any operator linear in heavy fermions will vanish when 
sandwiched between heavy sector vacuum. 

Combining the reduction processes of the bosonic and the fermionic sectors, it is 
clear that the connection term can be equally reduced to 

'-dLgMNZ^'^p'^tp'' ^ '-dsGaisy^tp'r ■ 

Now we can revert from the Hamiltonian to the Lagrangian, after putting all the 
heavy modes to their ground states, and arrive at the following reduced Lagrangian, 

•^for index only 
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= i^G.^ry'' + 'i^VG^Av + r^,y>^) - A,r + 'i^j'^^rr , (b.s) 

where we introduced the notation, also used in the main text, J-" = A^* (rfW) and its 
gauge field A. We already defined G as the appropriate block of g, but now valued at 
M-n- In other words, G = Ai^ig)- Remaining fermions live in the co-tangent bundle 
of Ain, so the resulting Lagrangian is A/" = 1 non-linear sigma model on 2(n — 1)- 
dimensional manifold A^n, coupled to an Abelian gauge field W. Supercharge of this 
dynamics is a Dirac operator on Al„ coupled to Abelian gauge field A, and therefore 
the index is, under the canonical choice of the chirality operator, 

Ch{J^)A{Mn) ■ 

Mn 
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